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ABSTRACT 


The bituminous concrete plant was studied and various 
methods of modeling the plant were reviewed and evaluated. 
The Monte Carlo simulation technique was selected and applied. 
A digital computer program was developed which models the ag- 
gregate gradation produced by the plant. 

The simulation first randomly selects the operating 
characteristics of a typical plant. The aggregate production 
is then randomly generated for the output of each of the sev- 
eral hot bins. The true gradation for each batch fe recorded 
together with the apparent gradation of periodic batches as 
obtained from the sieving of small samples. The apparent 
gradation contains errors due to sampling, splitting and 
sieving. 

The model can be used to apply different testing pro- 
cedures and specifications and study their effect. Thus, 
the desirability of various quality assurance procedures 
can be evaluated and an optimum testing strategy developed. 
The techniques employed and many of the building blocks de- 
veloped are applicable to other bulk material processes such 
as the portland cement concrete plant. 

The simulation described is a conceptual model based 
upon limited field data in New York. The State of New York 
has additional work in progress that may be used to further 
refine and calibrate the model. The program, coded in the 


FORTRAN language, is described in detail. 
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INTRODUCTION 


BACKGROUND 
General 

Bituminous concrete plants are used to produce mate- 
rial for highway pavements and many other surfaces such as 
shoulders, driveways, parking lots, airfields, aprons, dams 
and irrigation canals. In the highway field bituminous con- 
crete is one of two competing materials for high-type pave- 
ments and various asphalt mixes are used almost universally 
for the construction of intermediate and lower type roads 
and for the resurfacing, maintenance and patching of all 
types of paved roads and streets. 

Bituminous concrete plants are spotted throughout the 
countryside, often located adjacent to a quarry or fravel 
bank as a source of aggregate. As transportation is a sig- 
nificant factor of the total cost of these materials, a 
large number of plants have been built to supply local needs 
at a reasonable travel distance. 

The major components for these plants are manufactured 
by several firms, but the manner in which they are assembled 
and installed vary G7a819? Various Gert of the plants are 
constructed locally and adapted to the unique characteristics 
of the site and the requirements of a particular operation. 
Thus, even though all snente are basically similar, a wide 
variation exists between the capacities and their physical 


and operational characteristics. 


The plants may be privately owned or may be owned by 
a local governmental unit. In New York State some 120 plants 
are producing material used in State work. The State has a 
strong interest in controlling the operation of these plants 
to assure that a good product is produced. 

The final objective in the control of a bituminous 
concrete plant is to obtain the best pavement that is con- 
sistent with economic constraints and that will meet the en- 
vironmental and traffic loads to which it will be subjected. 
A good product is one that can be produced efficiently and 
economically and will perform well during a long, useful 
life. In addition, a good product should be uniform and con- 
sistent to permit periodic tests to be meaningful and to mini- 
mize adjustments in plant operation and in the placement and 
compaction of the material. 

The adequacy of the material in place is dependent 
upon several factors. The material properties are a function 
of the physical and chemical properties of the aggregate and 
bitumin and the surface characteristics of the aggregate. 
Plant operation variables effect the dryness of the aggregate, 
the gradation of the aggregate blend and the proportioning of 
the aggregate and bitumin. The properties of the mix are also 
dependent upon mixing time and temperature. The final product 


% 


is further altered by handling and construction practices (19) 


*Numbers in parentheses refer to entries in the List of 
References 


Description ofthe |Plant 


In order to understand the problems associated with 
modeling the plant, it is advisable at this point to briefly 
describe its makeup and operation (17). The components of a 
typical asphalt plant are illustrated in Figure 1. The type 
of plant shown is the batch mix plant which is presently the 
more common type in use. The other type is called a contin- 
uous plant, in which the weigh hopper is replaced by feed 
mechanisms which proportion the aggregate. The aggregate mix 
is then fed continuously into a pugmill where the asphalt is 
injected and in which the mixing takes place. In this study, 
the batch type plant was used to illustrate the process, but 
the model developed is applicable to both types of plants. 

The production process for the batch mix plant may be 
described as follows: 

(1) The aggregate is transported from separate stock- 
piles by truck or bucket loader and is placed into 
individual cold bins. As the material in the stock- 
piles is sorted by size, the aggregate in the cold 
bins is similarily sorted. 

(2) The aggregate is volumetrically dispensed from each 
cold bin onto a common feed belt. At this point, 
the aggregate sizes become intermixed. 

(3) The cold aggregate blend is then fed into a dryer. 

(4) The dried Mbomecate blend is transported to the 
gradation control unit at the top of the plant 


tower. 
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(5) Vibrating screens separate the aggregate into 
size ranges. 

(6) The heated and sorted aggregate is held in the 
hot bins. Plants generally use three and some- 
times four hot bins. 

(7) The aggregate blend is proportioned by with- 
drawing the various size ranges from the hot 
bins. Sampling apertures are located below each 
hot bin. The material from each bin may be sam- 
pled by holding a sample tray or a shovel in the 
stream as the aggregate drops from each of the 
bins. 

(8) The aggregate drops into the weigh hopper where 
it is proportioned by weight to conform to the 
prescribed plant mix ratio. 

(9) The blended aggregate then drops into the pug- 
mill where asphalt is added and the material is 
thoroughly mixed. 

(10) The mix is finally dropped into a truck for trans- 

portation to the job site. 

Drying of the aggregate is required to properly coat 
the aggregate with bitumin, thereby Aeassccine the desired 
physical properties of the asphalt mix. Secondly, the dried 
aggregate can be more precisely sorted. If screening of the 
cold or wet aggregate eee to be attempted, many of the smaller 
particles would adhere to the large particles and be carried 


over to the larger sized bins. Moreover, the wet material 


would tend to clog the openings in the smaller-sized screens. 

It may be seen that the function of the feeds from the 
cold bins is to closely approximate the desired aggregate 
blend to assure that the proper ratio of sizes is fed to the 
tower. If an insufficient quantity of one size is fed through 
the dryer, one of the hot bins may become empty. Conversely, 
if too much of one size is supplied to the gradation control 
unit, one of the screens may become overloaded, which results 
in finer material being carried over into a coarser bin. In 
either event, the end result disrupts the planned aggregate 
gradation (9). 

Gradation Control 

Various user groups have established limits within 
which the distribution of aggregate particle sizes must lie. 
These limits are based in part on the theoretical basis of 
providing a dense well graded blend in which the smaller 
sizes nest between the voids left by the larger sizes. The 
limits are also based on experience gained by correlating 
the size distributions used and the resulting performance of 
the pavements. 

Table 1 lists the New York State sieve series used in 
determining all gradations. Figure 2 illustrates the grada- 
tion specification currently used by the New York State De- 
partment of Transportation for top coarse mixes, type 1A (12). 
It shows the band within which the distribution of particle 
sizes must lie. In using this specification, a contractor 


selects a job mix formula (JMF) which passes through this band. 


NEW YORK STATE 
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Sieve 
Size 


1/2: inch 
1/4 inch 
1/8 inch 
20 mesh 
4Q mesh 
80 mesh 
200 mesh 
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FIGURE 2. BITUMINOUS CONCRETE GRADATION SPECIFICATION 


The formula chosen is a function of the particular aggregate 
used and the characteristics of the plant. Having selected 
this job mix formula, he is then held to producing an aggre- 
gate blend that does not vary from it more than a prescribed 
amount. These tolerances were selected to exclude the blend, 
if the percent of material retained on any sieve varied more 
than two standard deviations as determined by a comprehensive 
field survey and study previously conducted by the State (5). 
The contractor establishes a plant mix ratio (PMR), 
which is the percentage of aggregate that is to be taken from 
each hot bin to form the final mix. It is based on the ex- 
pected gradation in each bin and the desired job mix formula. 
Modern equipment, especially for those plants with automatic 
controls, permits close adherence to the specified plant mix 
ratio. Errors are likely, however, due to the material in 
the hot bins varying from the desired distribution of sizes 
and due to the segregation of sizes within the bin itself. 
The actual distribution of sizes in the resultant mix 
can be estimated by taking a sample of the final mix. The 
bitumin must be removed through an extraction process before 
the aggregate is graded. The aggregate is then passed through 
a series of sieves and the percent passing or retained on each 
sieve is measured by weighing. An extraction of the bitumin 
requires an hour and a sieve analysis can be done in about 
15 minutes. Batches, on the other hand, are produced at a 


rate of about one per minute. 
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The standard testing procedure for an asphalt plant 
is to sample at each hot bin as the material is dropped into 
the weigh hopper. These samples, one for each hot bin, are 
then sieved to find their gradations. The results are then 
combined mathematically using the plant mix ratio to esti- 
mate the overall ered aged cri of the mix. This is called a 
complete hot bin analysis and requires about 45 minutes. 
Even if this procedure were done for each batch, by the time 
the test was performed and analyzed, the batch would most 
likely be in place on the road surface. Because typically 
such complete sieving tests are only performed about once an 
hour, this procedure merely provides historic data. 

In an attempt to reduce the testing effort required, 
New York State has instituted a simpler test known as the 
hot bin uniformity test (13). The sampling procedure is that 
already described, but the sieving procedure is simplified. 
Each hot bin has primary size, defined as that size range in 
which most of the particles fall. For example, for the number 
1 bin the primary size may be material that has passed the 1/2 
inch sieve and retained on the 1/4 inch sieve, whereas the 
number 2 bin has material between 1/4 and 1/8 inches as the 
primary size. The uniformity test determines only the percent 
of the primary-sized material for each bin. While the test is 
somewhat simplified in that fewer sieves must be weighed and 
fewer computations are required, the time required is still on 


the order of 20 minutes. 
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PURPOSE AND SCOPE 


Research Objectives 


The long term objective of this research is to improve 
the quality assurance of asphaltic concrete. The work is 
limited at this stage to providing a desired blend of aggre- 
gate sizes. It should be noted that this alone does not nec- 
essarily assure a good pavement. The problem of optimal mix 
design is a most important but clearly separate problem. We 
are concerned here only that, given a specified aggregate 
blend, how can we best assure that it is produced within 
acceptable limits. 

The specific objective of this individual project was 
to show the feasibility of modeling an asphalt plant and to 
demonstrate the applicability of the simulation method. It 
was not considered possible in this first project to produce 
a highly sophisticated and completely validated model. It 
was felt, however, that a preliminary model would show the 
usefulness of the simulation method and demonstrate the 
application of the model to problems such as the relationship 
between testing strategy and process control. 

As more comprehensive field data are collected and 
analyzed, the model can be further refined and calibrated. 
Such a model then could be extensively run under controlled 
laboratory conditions to test many control techniques and to 


develop an optimum control strategy. 


=i? = 


scope of the Study 

The purpose of the project was limited to the study 
of the aggregate blending process in the bituminous concrete 
plant. Specifically, the research tasks included: 

1. A review of the state-of-the-art. 

2. An analysis of existing data. 

3. The design of a conceptual model. 

4. The development of a computer program for the 

model. 


§. Showing the applicability of the model to achiev- 


ing the research objectives stated above. 


SOURCES OF DATA 


The New York State Engineering Research and Development 
Bureau has conducted extensive field surveys of bituminous 
concrete plants. Analysis of these data gives some insight 
into the process that occurs in the plant and is necessary 
in order to model the operation. Observations must also be 
used to select parameters and calibrate the model that is 
developed. 
1961-1964 Survey 

During the period 1961 through 1964, 55 plants pro- 
ducing top course mix were visited at which 868 hot bin and 
682 mix samples were obtained (5). Of these plants, 51 were 
batch type and 4 were continuous plants. At this time auto- 


matic control was just beginning to be instituted and at 16 
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plants batching was performed manually, 32 were semiauto- 
matic and 7 were automatic, including the 4 continuous 
plants. Batches were sampled about once each hour. 

The initiation of this testing program focused atten- 
tion on the problems and apparently resulted in decided im- 
provements, as the variations were markedly reduced after 
the first year. Therefore, the 1961 data were omitted from 
further analysis. Other data were then omitted for various 
reasons, such as the number of observations being too small 
to be significant. 

The remaining data that were analyzed for this project 
then consisted of 341 samples from 20 plants collected during 
the period 1962 through 1964. Each of these samples yielded 
the gradation for each hot bin as sampled when it dropped into 
the weigh hopper. The percent retained on each sieve was mea- 
sured and the means and standard deviations were computed. 

As these samples were taken at intervals of approximately 
an hour, the results give no information on the short-term 
variations for which the plant can be assumed to be operating 
steady-state with a relatively constant input and under stable 
conditions. Furthermore, only the output was observed, and no 
information was gained about the internal. occurrences within 
the plant. 

1967 Survey 

As part of this project, technical assistance was given 

to the State in setting up a pilot study that was conducted 


during the summer of 1967. This was a field survey in which 
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intensive sampling was performed specifically to gather the 
type of information needed to model the plant process. 

Samples were taken at the beginning and then the end 
of the process separated in time by about 10 minutes. As 
this is the time required for aggregate to move through the 
plant, the assumption can be made that the same material is 
observed at both ends. 

Pairs of samples were also taken at each end as close 
in time as physically possible. Thus, the populations from 
which the pair of samples were extracted can be treated as 
constant; and the differences observed between the pair of 
samples is indicative of the errors inherent in extracting 
a small sample from a large flow. 

The sampling procedure was as follows: 

1. Three samples were taken at the cold feeds, one 

for each bin. 

2. Immediately afterwards, a corresponding set of 

three samples were taken at the cold feeds. 

3. Ten minutes later, three samples were taken at 

the hot feeds. 

4, Immediately afterwards, a corresponding set of 

three samples were taken at the hot feeds. 
An entire set then consists of 12 samples. This work was 
performed at the continuous plant at Middeville, New York. 
Thirty sets of 12 samples each were obtained. 


1968 Survey 


Based on the experience gained during the 1967 survey 
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a large scale testing program was established for 1968. 
Detailed field observations were made at 46 bituminous con- 
crete plants around New York State. For 39 of these plants, 
sampling was performed on a comprehensive basis. Four plants 
were intensively sampled using the same procedure employed 
for the 1967 study. The processing of the resulting 15,000 
samples (drying of cold samples, splitting and sieving) has 
been a major undertaking requiring 1 1/2 years. 

This study represents the first extensive time series 
data ever acquired for a bituminous concrete plant. The anal- 
ysis of the results will provide the information needed to 
calibrate and refine a simulation model. This analysis was 
not available, however, for this present project. The model 
was based on the 1962-1964 and the 1967 surveys. Nevertheless, 
the model was constructed to facilitate refinement when the 


results of this extensive 1968 survey are completed. 
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DEVELOPMENT OF A MODEL 


THE CONCEPT OF MODELING 


A model is an approximation to reality. Usually the 
closer it emulates reality the more complex it becomes. Such 
complexity makes it more difficult to use as an aid in inter- 
preting or understanding the process it models. Hence, there 
is always a trade off between Simplicity and validity, with 
the criteria being to make the model as simple as possible 
and still perform the htaskprequired. (2.8). 

Physical Models 

There are several modeling techniques which might be 
applied to this problem. The first is the physical model 
which is a direct representation of the process. A full- 
scale model would be the asphalt plant itself. One could be 
used specifically for research purposes. To provide better 
control of the aggregate flow process the material could be 
recirculated; thereby providing uniformity of input and a 
steady-state situation. Some degradation and deterioration 
of the aggregate would occur under these conditions in which 
particles would be broken down into smaller sizes and very 
small particles would be lost as dust. 

The Barber-Greene Company conducted such full-scale 
tests to study the operational characteristics of its dryer 
(1). Such an approach, however, is costly and it is diffi- 


cult to control the many operational and environmental factors 
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which affect the process. 

The overriding difficulty with this approach is that 
one is faced with a process that is most difficult to ob- 
serve, in which the sampling must be done sporadically in 
very small quantities and for which the testing procedure is 
tedious and time-consuming. 

An alternative approach is the small-scale physical 
model, which may be brought into the laboratory where many 
conditions may be more precisely controlled. In the asphalt 
plant, however, we are dealing with an aggregate flow in 
which the particles already range down to extremely fine 
Sizes. Miniaturization offers no real advantage and would 
be extremely costly. 

Analogous Models 

A second method would be to develop an analogous model. 
Stephens investigated the problem of aggregate segregation 
within a bin by utilizing colored ball bearings of a uniform 
size (16). While such an experiment might give some insight 
into the real problem, the obvious shortcoming is represent- 
ing the full range of sizes found in an aggregate flow. The 
provision of balls of varying size puts us right back into a 
full sized physical model. 

A different approach would be to model the aggregate 
flow by representing the individual particles by numbers in 
an array in a digital computer. Such a procedure has been 
employed successfully for the simulation of vehicular traffic 


flow. In simple form, a vehicle has been represented by a 
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"1" and a space by a "0" in a binary word (4). A more 
sophisticated procedure used an entire computer word to 
represent a vehicle by coding its location and character- 
istics into theoword (7)% 

Such a digitized analogous model was investigated 
for this proreeis wht ee brute force technique in which 
individual particles or particle groups are processed by 
the computer. Individual components such as the screens 
and the bins are represented by building blocks through 
which the aggregate flow is moved. For example, Figure 3 
is an array which simulates the flow in a bin. The flow is 
treated two dimensionally by considering only one plane 
through the bin. Although the concept could readily be ex- 
tended to three dimensions, the magnitude of the problem 
would be tremendously expanded. 

In Figure 3, seven different particle sizes are rep- 
resented by the digits "1" thru "7". The program removes 
digits from the bottom and then scans the bin, row by row, 
from the bottom up. Digits are allowed to drop down into 
voids created in the row below. Then new digits are entered 
from the screens above. Rules based on probability theory 
are used to determine how individual digits flow thru the 
bin. For example, a digit may drop directly downwards or on 
either diagonal. A digit near the center of the flow may have 
a higher probability of dropping than one near the wall (shown 
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This procedure is a direct analogy for flows in which 
all the particles are the same size. It must be modified to 
account for the aggregate problem. Basically each digit must 
represent a unit volume of a particular size. The rules for 
processing the flow then give the opportunity for the smaller 
particles to sift through the larger ones. It should be noted, 
however, that some 125,000 particles passing the 200 sieve 
could fit into the same box as one 1/4" particle. The lim- 
itations and inefficiencies of this model are readily appar- 
ent. 

Functional Models 

A functional model describes the basic interrelation- 
ships among the variables and parameters. In essence, the 
process is reduced to a mathematical equation. To accomplish 
this, much knowledge is needed about the entire process. As 
previously mentioned, little is known about the internal flow 
within the asphalt plant and the number of variables is very 
large. This procedure was considered to be beyond the capa- 
bility of the state-of-the-art. 

Prediction Models 

A more relaxed form of the functional model is some- 
times referred to as a prediction model - one that permits 
inferences to be made without complete knowledge of the inner 
workings of the system. For management purposes such models 
often suffice (3); 

One approach is to develop a prediction equation through 


the use of regression analysis. The output of the plant would 
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Be observed as. well as the input and various plant process 

and environmental factors. Factors to be recorded would in- 
clude: dryer type, size, angle charging rate, burner nozzle 
and air input; screen area, wear, speed and loading rate; 

feed mechanisms, type and settings; aggregate moisture con- 
tent, type, shape and gradation; plant mix ratio, number of 
batches since plant mix ratio was changed and number of batches 
Since beginning of day; temperature, humidity and time of day. 
Through regression techniques the output would then be cor- 
related with the significant variables and a prediction 
edUation established i(20): 

For the asphalt plant the situation is considerably 
complicated by the fact that the output cannot be described 
by a single dependent variable. The raw data consists of 
percent of material retained on each of eight sieves for 
three or possibly four bins. These may be combined by using 
the plant mix ratio to calculate the percent retained as if 
the sample were taken from the pug mill. We still have eight 
dependent variables which are in themselves interdependent in 
Siete they must all total 100 percent. That is, 1f one percent 
retained increases, one or more of the remaining ones must 
show a corresponding decrease. 

One solution to the problem of multiple output vari- 
ables is to convert them to a single index to use as the de- 
pendent variable. One index that is often used is the sum of 
the squares of the deviations. For each sieve, the deviation 


between the observed percent retained and the job mix formula 
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is determined and then squared. Figure 4 illustrates these 
deviations. The sum of these quantities for all sieves is 
the single index of how well the output compared with the 
desired gradation. The smaller the index the better the fit. 

The problem with this index is that deviations for 
some sizes may be more sel A Soba: than deviations for other 
Sizes. Furthermore, negative deviations may be more critical 
than positive deviations. A system of weighted deviations 
could accommodate these criteria. Further work would be 
needed to establish these weights. 

The major problem with the single index approach, in 
addition to it being a significant change from standard prac- 
tice, is that it averages deviations over the range of sieve 
sizes. The accepted procedure is the more comprehensive re- 
quirement that specifies that gradations must fall within 
an envelope. Quality assurance specifications frequently 
set. different.toleranees for Gach™percent retained. ~ However, 
for the purposes of comparison the somewhat crude Single index 
for the goodness of a gradation curve has considerable appeal 
and offers definite advantages. 

The development of separate prediction models for each 
of the output variables as functions of various plant and en- 
vironmental factors, would also beua possibility. If appre- 
priate data were available, regression techniques should be 
helpful in the construction and evaluation of such models. 
Such data are not available currently and might be most dif- 


ficult to obtain at any time because of the difficulty (i) 
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to measure some factors, (11) to quantify certain influen- 
tial factors, (iii) to vary certain factors over a reason- 
able range in a production facility, and (iv) to observe 
variability in the aggregate characteristics and the envi- 
ronmental factors over a\short sampling period. The prob- 
lems mentioned could be resolved and an appropriate sampling 
plan designed, but the operational problems of carrying out 
such a plan would be quite large. 

Monte Carlo Simulation 

The modeling of the asphalt plant is inherently a prob- 
ability problem. The flow, mixing, drying, screening, segre- 
gation, and feeding of the aggregate will cause the gradation 
to vary from batch to batch. Moreover, the sampling and 
testing procedures introduce random errors. These variations 
may be handled as fa prageapility proplem (2), 

The modeling, however, cannot be handled rigorously 
and analytically by existing probability theory. Insuffi- 
cient information about the process exists. Such factors as 
the interdependence of the dependent variables and drift, 
which may or may not be accompanied by changes in control, 
preclude the direct mathematical approach. 

Monte Carlo techniques use a series of random numbers 
from analytical or impirical distributions to provide the 
random nature of the aggregate flow and the testing. A series 
of rules are followed to generate the aggregate distribution. 
A testing strategy may be specified and the test results dis- 


played. The amount of computation is so voluminous, however, 
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that this method would have been unthinkable prior to the 
advent of the high-speed electronic computer. 

A special advantage of this technique is that all 
variables may be precisely controlled and the entire pro- 
duction can be observed and analyzed. This work can be 
performed very rapidly in the computer, whereas comparable 
field observations would take years. In the case of the 
asphalt plant, the aggregate gradation is generated by the 
computer for each batch and its true gradation is continu- 
ously known. One can also measure the gradation through a 
sampling and testing procedure which contains inherent errors. 
Thus, the test results for selected batches are readily com- 
pared with the generated "true" gradations for all batches. 

A second important advantage over "real world" methods 
is the ability to reset the random number series for addi- 
tional runs in which different testing strategies are tried. 
By repeating the series the same identical aggregate flow is 
repeated. Therefore, the differences in results are inde- 
pendent of any differences in the observed flow. Of course, 
data from the real world are needed to validate models devel- 


oped for the suggested computer experimentation. 


FORMULATION OF THE MODEL 


Investigation of the problem of model formulation lead 
rapidly to the fact that very little was known about the in- 


ternal workings of some of the components of the asphalt plant 
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and about the interactions between these component parts. 
Little data have been collected at the cold end of the 
process because of the necessity to oven dry the sample be- 
fore a reasonably correct sieve analysis can be made. The 
process up to the hot bins is inaccessible and hostile. The 
high heat, dust, a ea ae machinery and elevated location 
make any inspection dangerous, even if the facility could be 
modified to provide for sampling. 

The only location at which sampling is provided for 
is at the drop below the hot bins, and even this is an awk- 
ward, and, difficult location, at. many plants:... This is, the, one 
point in the process for which considerable data are avail- 
able. Material can be sampled from the truck into which the 
pug mill discharges, but testing of such samples for the 
aggregate gradation is complicated by the fact that the bi- 
tumin must first be extracted. 

It. became clean then that the output of the. plant 
could be conveniently thought of as the drop from the hot 
bins. In accepting this one makes the assumption that the 
proportion of the output from the various hot bins is accu- 
rately blended in accordance with the prescribed hot mix 
formula. As typically this blending is done by weight with 
calibrated and automatically controlled equipment, the assump- 
tion of insignificant error beyond the sampled point is rea- 
sonable. 

The technique employed in formulating the model was 


then merely to simulate the output of the plant. This output 


pa 


was considered to be observed at the drop below the several 
hot bins. This procedure avoided many pitfalls that would 
have resulted in attempting to simulate the entire process. 
Furthermore, the plant simulated would be a typical New York 
state plant rather than a specific plant. This general simu- 
lation results from the feet that the modeling and calibra- 
tion is done with state-wide data collected from many dif- 
ferent plants. 

In the typical simulation, time is one of the funda- 
mental variables. The program processes the flow for an 
increment of time and then the procedure is repeated for the 
next time increment. The asphalt batch plant by its very 
nature has incremented time with the series of batches, and 
the continuous plant can be considered to be a stream of 
batch-sized flow elements. 

In the adopted procedure, time is treated indirectly 
as an implied variable with the individual batches handled 
consecutively. While most simulations require several cycles 
to load the system and reach a steady-state situation, this 


model is not dependent on any such initializing period. 


EVALUATION OF THE PROCESS 


Types ofurlrests 


Several testing methods may be employed to evaluate 
the aggregate gradations produced by the asphalt plant and to 
assure quality control (18). Three. such tests are described 


to demonstrate the testing procedure. 
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The tolerance test checks the percent passing each 
sieve for the combined mix sample representing the entire 
batch. This percent passing must lie within the range of 
the values deliniated by the job mix formula, plus or minus 
an.associated: tolerance limit. If, a batch has, a.percent 
passing outside this range for any sieve, the batch is 
deemed to be out of specification. 

The present approach used by New York State is to 
establish the tolerance limits so that five percent of the 
material evaluated by each sieve will fall outside the ac- 
ceptable range. As the area under the normal curve between 
-2 and +2 standard deviations is 95 percent of the total 
abe a » cthewtolerance i mags) Meréueat at +t) etandanrd devi ations, 

The uniformity test examines only the percent retained 
on the primary size sieve for a bin. It is a shortcut test 
as compared with a complete sieve analysis. The primary- 
size percent retained must not be less than the specified 
uniformity test limit. The value currently used by New York 
State is 70 percent for both Bins 1 and 2. If the primary- 
size percent retained is less than this limit for either bin, 
the batch is out of specification. 

The uniformity difference test monitors variation in 
primary size between the batches tested. The test compares 
the percent retained on the primary-size sieve with the pri- 
mary size of the previous batch tested. If the absolute 
value of the difference is greater than the specified uni- 


formity difference limit, the batch is out of specification. 
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The number of such checks is always one less than the number 
of batches tested, as the first batch tested cannot be eval- 
uated. New York's uniformity difference limit is presently 
12: percent, 

Testing Errors 

There are several accidental errors that are inherent 
in testing procedures used to estimate the gradation of the 
aggregate. The three primary sources of error result from 
Pemp lin Sop atti tard ea ey ine operations? (6 ).°¥ The? com 
bined error from these sources is herein called the "testing 
error". 

The sampling error occurs because a very small sample 
on the order of a gallon is used to estimate the character- 
istics of an entire batch. The magnitude of the sampling 
error may be estimated ine comparing samples taken very close 
together time-wise, so that the bin gradations may be assumed 
to be in a steady-state condition. 

Once a sample is obtained it is commonly put through a 
splatter toPobtain® the ‘amount? ofematerial: which? canbe con- 
veniently sieved. i'Some error is introduced as the splitter 
does not divide the sample into two parts with identical gra- 
dations. When four quarters of a split sample which was split 
twice are sieved, however, the comparison of the results can 
be used to estimate the splitting error. 

A further error is introduced during the sieve analysis, 
because sample sieving is never completely reproducible even 


when the same sample is resieved. This error can be determined 
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by the sieving, recombining and resieving the identical 
sample. 

The testing error represents the difference between 
the true gradation of a batch and the estimate of that gra- 


dation as obtained by sieving a sample taken from the batch. 
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Testing Strategy 


A testing strategy is the sequence in which the var- 
ious tests are applied and the decision rules that are uti- 
lized to evaluate the material produced by the plant. For 
exampile ,..one may elect sto’ sampler everyr tenth batch: and to 
employ three uniformity tests followed by one tolerance test. 
If, however, a batch fails the uniformity test, a special 
sample will be taken and a tolerance test applied before 
deciding to accept or reject the product. Such a strategy 
involves a specific testing effort. A goal that can be 
studied by use of the simulation model is the trade-off 
between effective control of the process and the testing 
effort applied. 

One criteria that may be used to evaluate a testing 
strategy is the number of batches that are outside of speci- 
fication. The basic criteria of control is to reduce the 
variation of the product; therefore, the testing strategy 
should reject the tails of the production distribution while 
accepting the bulk of the batches which lie close to the mean. 

Due to sampling inadequacies and the testing error, 
some material may be rejected that is in reality acceptable 


(Type I Error) and, conversely, other material may be accepted 
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that should be rejected (Type II Error). It is desirable to 
keep these so called producer's and buyer's risks within 


acceptable bounds (10). 


ave 
DESCRIPTION OF THE; MODEL 


THE PROGRAM 


The basic concept of the simulation is quite simple 
and easy to comprehend, if one does not get lost in the 
complex details of programming. In essence, a typical 
plant is selected and batch by batch the production for, 

Say, one day is generated. The true gradation is recorded 
(an item never known in the real world), as well as the 
apparent gradation as obtained from the sieving of a small 
sample from each batch. The results of various tests are 
then determined for batches at selected intervals, and the 
action taken in accepting or rejecting batches is noted. 

The effect of these sporadic tests of the apparent grada- 
tions are then compared with the proper decisions that would 
have occurred from knowing the true gradation of all batches. 
Production 

A step-by-step description of the program is presented 
in ‘the simplifiedsfilow\ehabtiy shown as tipure o. The first 
step is to initialize the program and read in the input vari- 
ables..., Theseyvankabbes;itustedivini Table 2, describe the 
plant to be simulated, give operating instruction to the pro- 
gram, calibrate the process based on field observations and 
set forth the specifications and the testing procedure to be 
followed. The specifications are described in detail in 


Appendix D. 
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Input and Initialization 









Generate % retained for 
each sieve and each bin 


Determine testing errors 


yes 


ie 2 Introduce drift 


Calculate % passing all bins 
Combine using PMR and 
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Compute JMF as mean of production 


JMF acceptable ? 


yes 
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and perform statistical analyses 


Display output 
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FIGURE 5. BASIC FLOW CHART FOR THE SIMULATION MODEL 
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INPUT INFORMATION 


Plant Description Variables 
Number of bins ) 


Number of sieves 
Planti mix ratwe 


Primary size for each bin 


Program Operating Instructions 


Number of plants to be simulated 
Number of batches per plant 
Frequency of sampling 

Amount of output to be displayed 
Actuator fon. dri tt.option 


Random number seeds 


Calibration Data 


Mean and standard deviation of percent retained on 
each sieve for each bin 


Mean and standard deviation of testing errors 
(adjustments) for each sieve for each bin 


Testing Procedure and Specifications 


Upper and lower limits of acceptable job mix formula 
for each sieve size 


Tolerance test limits for each sieve size 
Uniformity test limit 


Uniformity difference test limit 
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The next step is to generate the percent retained on 
each sieve for each bin. This is accomplished by generating 
a random fraction and using it to enter a cumulative frequency 
plot’ for each sieve.:*Figure:6 illustrates this procedure, 
whereby the random fraction 0.56 is used to obtain a 77 per- 
cent retained on the 1/8 inch sieve of Bin 2. The curve is 
stored as a series of points and linear interpolation is 
used to obtain the actual value. At the same time the stan- 
dard deviation of the percent retained and the mean and 
standard deviation of the testing error are obtained. These 
three values are stored as a function of the percent retained 
and are selected by the same random fraction. A more de- 
tailed description of the procedure for storing these vari- 
poleseds tshown Gin Pulgure Ss "or AppendisoD.. 

The adjustment to be applied for the testing error is 
the mean adjustment plus a randomly selected standard devi- 
ation for the adjustment. The adjusted percent retained is 
then the true percent retained plus the adjustment. 

The percent retained for all sieves must total 100 
percent, because none of the sample is lost. However, as 
the percent retained values are selected randomly for each 
Sieve, there is no built-in balance to force them to total 
the necessary 100 percent. Hence, these numbers are stan- 
dardized in effect by multiplying each value by the ratio 
of 100 to the sum of the individual values. It is possible 
under this mathematical procedure for one of the factored 


values to become negative -- a physical impossibility. If 
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this occurs, the value is set to zero and the remaining val- 
ues are factored to total 100 percent. The question then 
arises as to whether the resulting standardized values have 
distributions that are similar to the original observations. 
Appendix B demonstrates that such an assumption is reasonable. 
The percent retained on each sieve is then combined 


for all bins and the percent passing is computed for the com- 


bined sample. The plant mix ratio gives the proportion of 
the final mix taken from each bin. Table 3 illustrates this 
calculation. 


The above procedure is then repeated for each batch 
until the specified amount of production is completed. The 
job mix formula is computed as the mean of the percent passing 
for each sieve for the entire production. This is a reason- 
able assumption as the job mix formula is typically based on 
the past gradation records of the plant. Field observations 
have shown that the job mix formula is a close approximation 
of the average gradation when a plant is operating properly 
C5). 

The individual values of the percent passing each 
Sleve as calculated for the job mix formula are checked 
against the State's band of acceptable values. If any sieve 
value is found to lie outside the band, the job mix formula 
is not acceptable and would not be approved. Therefore, the 
plant is not realistic and the simulation is terminated. As 
the band is quite large, this circumstance would rarely occur. 


An adjustment of the plant mix ratio corrects the problem, if 
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it should develop. 
Testing and Evaluation 

Dhewvariousictestscare then applied rte;the plant's ipro- 
duction. © Table th clists :thenoutput sinformatron which sis; idis- 
played for each plant that is simulated. 

Tolerance Test. The standard deviation of the percent 
passing is calculated for each sieve when the job mix formula 
1s Foundlevit ms multiplied thy \2utouform the icakeubated two- 
Sigma tolerance limit. These limits, however, are not used 
in the actual tolerance check, but rather the imposed toler- 
ance limits which were specified in the input are used for 
this purpose. The calculated two-sigma limits are displayed 
along with the imposed tolerances for comparison purposes. 

Each batch has its percent passing for each sieve size 
checked to determine whether it falls within the imposed two- 
sigma tolerance limits from the job mix formula. If a batch 
has any percent passing outside the limits, the entire batch 
is declared outside of specification. The tolerance test is 
the base on which the two other tests (uniformity and uniformity 
difference) are evaluated. Whatever the tolerance test labels 
a batch is assumed to be the true situation. 

Uniformity Test. The uniformity test utilizes the fi- 
nal percent retained on the designated primary size sieve for 
those bins having a specified primary size. Each percent re- 
tained is compared to the uniformity limit given in the in- 
put. If the value falls below the specified minimum, a reject 


is indicated. 
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OUTPUT INFORMATION 


Input Information 


(See Table 2 on page 34) 
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Job Mix Formula 


JMF calculated as the mean for the plant's production 


For-each sieve size the 


calculated and imposed toler- 


ance, upper and lower limits are displayed 


Tolerance Test 


Table showing the batch 
rejection 


Histogram of the number 
occurred 


Histogram of the number 
batches 


Uniformity Test 


Table showing the batch 
rejection 


Histogram of the number 
occurred 


number and sieve sizeuroreeacan 
of sieves on which a rejection 


of times each sieve rejected 


number and bin number for each 


of bins in which a rejection 


Histogram of the total number of rejects by each bin 


Uniform tyaDaetenence lest 


Batch numbers and bin numbers for each rejection 


Histogram of the number 
occurred 


of bins in which 4) relecuua: 


Histogram of the total number of rejects by each bin 


Summary of Results for all Three Tests 


For each test the number and percent of batches is 
given for each of the following: 


(1) Total acceptances 
(2) Total rejections 


(3) Acceptances that should have been accepted 
(4) Rejections that should have been rejected 
(5) Acceptances that should have been rejected 
(6) Rejections that should have been accepted 
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Uniformity “Difference Test. ~The uniformity differ- 
ence test utilizes the final percent retained on the primary 
Size sieve from batch to batch. Each percent retained is 
compared with the corresponding value for the previously 
tested batch. If the absolute value of the difference is 
greater than the specified input uniformity difference limit, 
a reject is indicated. 

Table of Results. Following the application of the 
various tests, the testing limits and strategies are evaluated. 
The number of tolerance, uniformity and uniformity difference 
test acceptances and rejections are tabulated and the percent- 
ages of acceptances and rejections by each test are calculated. 
At this point the tolerance limit for each sieve can be de- 
termined for the plant to yield the desired percentage of 
out-of-specification batches that is deemed acceptable. 

The uniformity and uniformity difference limits can be 
calibrated not only according to the percent of samples out- 
Side of specification, but also so that the Type I and Type 
II errors are at an acceptable level. At present there are 
no state guidelines on the acceptable magnitude of Type I 
and II errors, however, 10 percent risks for both the pro- 
ducer's and buyer's risks seem reasonable. 

Random Variable Generation 

A series of uniformly distributed random numbers is 
generated by the multiplicative congruential method in which 
each succeeding number is a mathematical function of the pre- 


ceding number. 
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To,illustrate let yNowbe .the .anttidgl spandommnumber 
or the seed of the number series. Subsequent numbers are 


then generated according to the formula 
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where e is a suitable multiplier, b is the number base of 
the computer and n is the number of digits in the word size 
of the computer. The resulting number is then divided by the 
largest number that can be generated to give a decimal frac- 
tion which is uniformly distributed over the interval be- 
tween zero and unity. 

The series may be reset at any time to repeat itself 
by setting N. qyback tome « Likewise, several independent 
random number series may be maintained by using different 
seeds. 

Random variables are atee required that are not uni- 
formly distributed, but follow some other prescribed distri- 
bution. This is accomplished by utilizing the fact that any 
random variable, no matter what its distributional form, has 
its cumulative distribution function uniformly or rectangu- 
larly distributed over the interval (0, 1); Ldenoted R (0,1)]. 
Therefore; by using a random fraction from R (0, 1) one can 
obtain the corresponding value for the argument x of the 
cumulative distribution function F (x) for the random vari- 
able — | thatais 0 xe F(x), the inverse of the cumulative 
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This sprocedure+is demonstrated in Figure 7.. Here, a 
uniformly distributed random fraction between 0 and 1 is used 
to obtain a normally distributed random variate between -3.4 
and +3.4, with a mean of zero. The cumulative distribution 
curve 1S approximated by a series of straight line segments. 

The result in this example may be used as a randomly 
selected deviation (d) to be applied to a known mean (x) and 
known standard deviation (SG ).- That is, any particular 


value (x) is computed as 
3) BS Beet Gos, 


In the case of a normally distributed variable a mathe- 
matical solution could be He See in place of the described 
method. In many instances, however, the distributional form 
is only graphically known. For example, Figure 6 illustrated 
the selection of a percent-retained value which is based on a 
distribution observed in the field. 

This procedure was used in the study primarily for its 
versatility. That is, for the same random number generation 
routine, any form of F(x) can be inserted be it empirical, 
analytical or an approximation to an analytical expression. 
It was considered that such generality more than offset-the 
loss of efficiency as compared with routines tailor-made for 
a specific distributional form. 

Drift Option 
A drift routine is contained in the program and is 


available as an option called for by the input control cards. 
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The provision of drift adds a realism to the simulation which 
could not be accommodated as easily in an analytical solution. 
It enables one to investigate a testing scheme's relative 
Sensi tinaiyein, detectang Grift® an the production process. 

Field observations indicated that gradations occa- 
Slonally tend to drift from the overall mean of plant pro- 
duction. This drift is primarily due to external forces 
affecting normal plant operations. Gradations drifted back 
to the mean either when corrective measures were instituted 
or when external pressures causing the drift subsided. 

Fluctuation in the rate of production is a major in- 
fluence on gradation uniformity. Production surges occur 
frequently during plant startup when many trucks are waiting 
to be loaded. As plant production approaches capacity, more 
aggregate must be separated by the vibrating screens and screen 
efficiency drops with a resulting carry-over of materials. 
This carry-over may occur when excessive charging of the 
screens causes smaller material to be trapped and carried 
along with the larger stone or when screens become plugged. 

In either case the aggregate is not allowed to drop into the 
proper hot bin. 

Torn or worn screens also influence production unifor- 
mity by permitting oversized particles to drop into the wrong 
hot bin. Carry-over and faulty screens tend to reduce the 
percent-retained primary size in the hot bins. This reduction 
may be considered as a negative drift from the plant mean and 


is the more prevalent type of drift. 
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An increase in the percent-retained primary size, 
termed a positive drift, is also possible.. It may occur 
when plant production slows to well below normal. Numerous 
other variables such as stockpile gradation and cold feed 
intermixing may also shift the mix gradation. 

The program only applies drift to the percent retained 
on the primary size sieve for bins which have such a primary 
size designated. All sieves for the bin are affected, how- 
ever, when the percent-retained values are adjusted to total 
100 percent. 

When the drift option is activated, the type of drift 
(negative or positive), the length of the drift in batches, 
the intensity of the drift and the frequency of drifting are 
selected as a random process. The drift applied complies 
with a parabolic function that reaches its maximum value in 
two-thirds of its duration and then returns to the mean. 
Subroutines 

The program was built as a series of building blocks 
to facilitate change and refinement. Several sub-programs 
were used for reoccurring portions of the program. These 
SUBROUTINES are described as follows: 

(1) Subroutine N@RMAL accepts a uniformly distributed 

random number between 0 and 1 and uses it to generate 

a normally distributed random value between -3.4 and 

+3.4 with a mean value of zero. A cumulative normal 

curve is specified as a series of points thru a DATA 


statement. The value is obtained by linear inter- 
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polation between these points, thus assuming that 

the curve is made up of a series of straight line 
segments. The resulting value is used as a randomly 
selected standard deviation for normally distributed 
data. Figure 7, previously presented, illustrated 
this procedure. 

(2) Subroutine GETVAL selects a mean percent re- 
tained and its associated standard deviation. These 
values are obtained by linear interpolation from a 
family of curves input to the main program for each 
bin and each sieve. 

(3) Subroutine GRAPH plots frequency histograms in 

a bar chart form. Each graph is automatically scaled 
to fit on a single page. Up to ten different fre- 
quencies can be accommodated. The arguments are the 
number of different frequencies, the individual values 
of each frequency, and the title of the graph. 

(4) Subroutine ZER®@UT sets the entries in any array 
to zero. It is used to zero the frequencies used in 
subroutine GRAPH. 

(5) Subroutine RANDU generates a uniformly distributed 


mandomerrpaction in the interval O°to 1. 


CODING 


The program was coded in the FORTRAN IV compiler lan- 


guage which is extensively used by engineers and scientists. 
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The use of this coding language facilitates communication 

and results in a product that may be understood, used and 

modified by other researchers. To aid in understanding 

the program,detailed information is given in the Appendices. 
Appendix ¢©-lists,and (defines;.the’ variables,~used. The 

actual program is listed\as Appendix F.. s-Thesfrequent) use of 

COMMENT statements in the listing serves to further describe 


the detailed aspects of the program. 
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FINDINGS 


CONCLUSIONS 


Within the state-of-the-art, the Monte Carlo simula- 
tion technique offers the most potential for modeling the 
bituminous concrete plant. Such a simulation can be ef- 
fectively performed by a digital is oh eons een, lit te 
sufficient knowledge exists to model the entire process 
within the plant, however, a model of the plant's output 
is both feasible and acceptable. 

The conceptual model developed is sufficient to dem- 
onstrate the usefulness of the program and its applicability 
to the study of quality assurance procedures. The effec- 
tiveness of the program can be further enhanced by additional 
refinement and field calibration. 

The simulation model offers significant advantages over 
conventional field studies in evaluating various testing 
strategies and selecting test specifications and limits. 
Production can be investigated in the laboratory at low cost 


under highly controlled and reproducible conditions. 


RECOMMENDATIONS 


The extensive time-series field data collected by the 
State of New York should be employed to further improve and 
calibrate the simulation model. 


The model should be used to evaluate various testing 
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procedures and to select an optimum testing strategy and 
realistic test Specifications and limits. 

Further development of the model should proceed 
together with a modest program for the collection and 
correlation of field data in order that the model be 
firmly based on reality. 

Refinements in the program are needed in handling 
drift and to accommodate the situation where changes in 
the plant control variables are made during a production 
ran, 

The techniques employed in this research and the pro- 
gram building blocks developed should be applied to the 
study of quality assurance for other bulk construction 


materials, such as portland cement concrete. 
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GLOSSARY 


The terms in this list are defined in the specific sense in 
which they are used in this report: 


Final Mix. The mix of aggregate sizes produced at the 
output end of the plant as the composite of the material 
dropped from the hot bins. 


Job Mix Formula (JMF). The target gradation expressed 
as the percent passing each sieve for which the producer aims 
in® nes pina] ams x. 


Percent Passing. The percentage of aggregate by weight 
passing a given sieve as computed for the final mix by using 
the plant mix ratio and the results of the sieve analyses for 
the individual bins. 


Percent Retained. The percentage of aggregate by weight 
retained on a given sieve as determined by the sieve analysis 
of a sample for an individual bin. 


Plant Mix Ratio (PMR). The percentage of the aggregate 
by weight which is dropped from each of the hot bins to make 
Np ene ianal mix. 


Testing Error. The difference in the true percent re- 
tained for the aggregate in a bin and the observed value 
which contains the composite errors resulting from the sam- 
pling, splitting, and sieving operations. 


Tolerance Test. A test in which the percent passing 
each sieve for the final mix sample is compared with the 
value delineated by the job mix formula plus and minus a 
tolerance limit for each sieve. 


Uniformity Difference Test. A test which compares the 


percent retained on the primary-size sieve for a bin against 
the value observed for the previous batch tested to ascertain 
whether the absolute value of the difference is less than the 
uniformity difference limit. 


UnlcoOrmury Test. A test which compares the percent re- 
tained on the primary-size sieve for a bin against a specified 
minimum value called the uniformity limit. 


APPENDIX A 


EVALUATION OF DISTRIBUTIONAL FORMS 


In describing empirical data it is always 
helpful to be able to model their distribution with an 
analytical expression., If such an expression is used, 
the evaluation of its efficacy is desirable. Many 
measures that provide such an evaluation are available. 
None of these are particularly powerful to discriminate 
among various competing models. However, thev do provide 
an objective procedure for casting doubt on a model. 

When such doubt does not seem warranted at the risk 
levels required, the conclusion is made that the proposed 
models seem to be worthy of consideration and use. Two 
of the more powerful testing procedures are described be- 
low. The first is a classical procedure in common use. 
The second is a natty developed procedure which, for 
testing for normality of the raw data or transformations 
of the raw data, is more powerful than the former more 
general procedure. It is the latter that has been used 
for evaluating distributional models in this study. 

All tests of distributional models involve 
measures of discrepancies between what is observed 
empirically and what would be expected if the postulated 
model was appropriate. For instance, if the data are 


subdivided into k intervals and O; denotes the number of 


data observed in the i-th interval, e. denotes the 
number of data expected in the i-th interval if the model 
is correct, then the statistic 
k 
ets (Co. - e.) 
; he - 
ie) =e ee 


Cus 
ae 


provides a measure of the goodness of fit of the model to 
the data. If none of the e; are too small, then this 
statistic is approximately distributed as a chi-square 

with degrees of freedom equal to k minus the number of 
independent uses made of the data to estimate the e,'s. 

For testing for a Gaussian or normal model where the two 
parameters of the distribution are estimated from the 

data, the number of degrees of freedom is k-3. The decision 
procedure in this instance is as follows. Compute x? and 
reject the hypothesis that the distributional model is 


mai: xe ; k-3 where the probability of observing 


correct if X 
a chi-square random variable with k-3 degrees of freedom 
that exceeds xe ; k-3 is a (o< a <1). Otherwise accept 
the hypothesis of the appropriateness of the distributional 
model on the grounds that the data provide insufficient 
evidence to support an alternative. 

The second test (used in this study) is an 
analysis of variance test for normality developed by 


Shapiro ‘and Wilk Ci@4), * THe tdéet etratistie for cthisa test , 


referred to as the W-statistic, is a function of the 


deviations of the ordered observations from a straight 
line obtained by fitting by the method of least squares 
the ordered observations to the expectations of these 
ordered observations (where the expectations are based 
on the assumption of a normal distributional model). Of 
eoUunse, “Liv tais desired to test for lognormality or 
normality of the square roots etc., such transformations 
would only need to be performed before employing the W-. 


statistic. If the n ordered observations are denoted by 


nh 


then, Wilk and Shapiro (14) introduce a statistic 


re 2 
Gore a. cae 
W= i=l 
n 
i=1 
wheresy" =! Yisn, and the a.'s are related to the expected 


1=1 
values of the ordered observations from the standard normal 
distribution. These a,'s are tabulated in (14) for n=2 
( 1) 50 and an approximation procedure for evaluating W 
is provided for larger n. The exact null distribution for 
Wis also tabulated (U44)i tom naa Glo. 

In (14) and (15) considerable attention is given 
to examining the power of W as compared with other 


Os 


"goodness of fit" criteria, the X° test included, and it 


is found to be more powerful against many alternative models. 


It was for this reason that the W-test was chosen for 
evaluating the goodness of fit of the proposed distribu- 
tional models in this study. 

The operational procedure for use of the W- 
statistic is as follows: Compute W and reject the 
hypothesis of normality if 


Oe ih * ages Stile 


where the probability of observing a value of W in excess 


os es - a; n based on a sample of size is 1 - a when the 


model of a normal distributional form is appropriate. 

The W-statistic was computed for the raw data, 
the natural logarithms and square roots of these data for 
each of three bins from 20 asphalt plants sampled by the 
New York State Department of Transportation in 1962-64. 

A Summary of the decisions available using a 10% level of 
Significance is reported in the following tables. Although 
the evidence is not overwhelming, the single distributional 
model (for the primary sizes) that appears to be most 


Suitable is that of the normal distribution. 


TABLE 5 
ANALYSIS OF DISTRIBUTIONAL FORM 
BASED ON BIN 1 


Using W Statistic at 10% Level of Significance 
Number of Plants = 20 


Primary Size = 1/4 
Distributional Sieve Size ; 
Form THe ree OM eee Ur ode eo 
No. of Plants Accepted by Each Distributional Form 
Normal 8 137 ids 12 e 
5q.- Root 7 16 ee y2 2 4 
Log , 2 17 127 5 1 


No. of Plants Accepted by Only One Distributional Form 


Normal 2 53 2 0 0 
Sd. KOeT i 0 0 0 
Log 0 0 0 0 L 


No. of Plants Accepted by All Three Distributional Forms 


All Three y. 16 da 5 0 Af 


Note: Sieves 1/2, #80 and #200 had very 
small percent retained values. 


TABLE 6 
ANALYSIS OF DISTRIBUTIONAL FORM 
BASED. ON BIN 2 


Using W Statistic at 10% Level of Significance 
Number of Plants = 20 


Primary Size = 1/8 
Distributional Sieve Size 
Form Weer iter tis 20 80 200 


No. of Plants Accepted by Each Distributional Form 


Normal 0 Ins 18 20 6 5 
5q. Root 0 14 1 Li, 5 
Log 0 15 14 45 e} 1 


No. of Plants Accepted by Only One Distributional Form 


Normal 0 ia 3 3 1. 
Sa. ROOT 0 0 0 0 
Log 0 3 0 0 


No. of Plants Accepted by All Three Distributional Forms 


Me, Lnree 0 7 14 15 3 = 





TABLE 7 
ANALYSIS OF DISTRIBUTIONAL’ FORM 
BASED ON BIN 3 


Using W Statistic at 10% Level of Significance 
Number of Plants = 20 
No Primary Size 


Distributional Sieve Size 
Form Be 920 Og aS 8 20 80 200 


oe 


No. of Plants Accepted by Each Distributional Form 


Normal 0 di Bal ZU 20 20 
Sq. Root 0 He Ig Bic 14 14 
Log 0 0 Rig) £5 14 mies) 


S 


No. of Plants Accepted by Only One Distributional Form 


Normal 0 0 2 4 mn 4 
Sq. Root 0 0 0 0 0 
Log 0 0 2 0 0 0 


‘S 


No. of Plants Accepted by All Three Distributional Forms 


All Three 0 0 7 14 13 v3 


Note: Sieves 1/2 and 1/4 had percent 
retained values close to zero. 


APPENDIX B 
STANDARDIZATION OF MEASUREMENTS 


In the simulation model used, a plant was 
represented essentially by taking observations from various 
"Dercent retained" distributions representing appropriate 
classes of material. When such numbers are put together, 
their sum is not likely to add to the desired 100. Hence, 
for purposes of this investigation, the summing to 100 
has been assured by standardizing these numbers by 
multiplying each observation by the ratio of 100 to the 
sum of the observations from each class. This standardiza- 
tion raises the question as to whether or not the 
standardized values have distributions that are similar 
to the original observations. The rather cumbersome algebra 
that follows in this appendix leads one to believe that 
to assume so will not lead one too far wrong. Although 
based only upon a comparison of moments, which is certainly 
not a proof, the comparison seemed adequate for the 


purpose. 
Ay 





Let yx), Pog Se Le Di. tole CEO reey en aller 
variable with mean m(x. ) and variance a Cx.) Tat. ae. 
we have k distributions that are independent of each 
other with potentially different means and variances, 
respectively. Consider the selection of a single observa- 
tion from each of these k distributions and consider the 


Statistic: 


CAO AD (1) 


Each of the y(x,)'s represent a percent of material re- 


tained by sieve of size denoted by the X'S. Because of 
sampling errors the sum k may not add up to 100. 
x y (x, ) 
i=l 


By weighting each y(x,) by 100/ yx), the resulting 
weighted values y' Cx.) will ees 100. The question 
then arises: Are the y' (x. ) representing the appropriate 
distributions? One thing that one would want would be 
for ECy' (x. )) = m(x.). This is approximately the case 


as 1S indicated below. 


By taking logarithms of both sides of (1) 


yielding 
k 
ln vy Css) = eee meee CX. In. 2 oviC xin 
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Suppose that o (x, ) is proportional to m(x., ) (which may 


be a very reasonable supposition). Then: 


2 Zz 2 k 
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(The above inequality results from the fact that 


V4 
m? (x, ) > = Cr mx, 1? 


1=1 

which follows from Schwartz Inequality.) 

Hence, the above discussion of the mean and variance for 

y' Cx.) gives one reassurance that the transformed y(x,)'s 

transformed to guarantee a sum of 100) have expectations 

and variances not drastically different from the original 
y(x,)'s. The above approximate solutions depend upon 

the following: 


i.e. the first two terms of the Taylor series expansion 
of the natural logarithm about the mean. 
Ee ey (x, ) = m(x,) 
Vin Sy Gx as 4 Viv (xed 
ay ~~ i Naat 
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From the above, it is difficult to appreciate the magnitude of 
these covariances. For certain combinations of Xs and ae, the 
covariance is going to be close to zero. However, for the x; 
values near the primary size, the covariance is going to be 


negative and not necessarily small. 


APPENDIX C 


DICTIONARY OF FORTRAN VARIABLES 


Program Symbol 
MAIN PROGRAM 


Meaning of symbol 


Arrays 

ACCREJ Calpe sacei dae by tolerance test end rejected 
byi uniformity. or difference) tese 

ALLRET Distribution of mean and standard deviation of 
percent retained and adjustments 

BACC Counter-accepted by both tests (uniformity or 
difference and tolerance) 

BINNAM Bin names 

BREJ Counter-rejected by both tests (uniformity or 
difference and tolerance) 

DRFMAX Maximum drift 

DUR Duration of -drity¢ 

FREQ Counter —- total tolerance rejects 

FREQ1 Counter -—- tolerance rejects by sieve 

FREQ2 Counter —- total uniformity or difference rejects 

FREQ3 Counter ~ uniformity) or difference rejects by bin 

IRN Random number seeds 

TSAR. Batch ab. which “drittistaces 

L@LIM The. State's,» lower job mMixeiLormiulasband 

MSTEPS Increment between samples (same as STEPS) 

NPTS Number of points in ALLRET 

PARRY Accept reject matrix 

PASS Calculated percent passing 

PMR Plant mie cae 

PRMBIN Primary size sieve designation 


Program Symbol 


REJ ACC 


RETAIN 
RN 
SNAME 
STATE 
STEPS 
TEMP 
TEST@L 
TOTTI 
TIT2 
TPT 3 
TIT4 
TITS 


T@LACC 


UNACC 


UPLIM 


ADJSIG 
ADJUST 
AJMEAN 
BATCH 


BIN 


* BLANK 


Meaning of Symbol 
Arrays 


Counter - rejected by tolerance test and 
accepted by uniformity or difference test 


Calculated percent retained values 
Generated random numbers 

Sieve names 

Imposed tolerance limits 

See MSTEPS 

Temporary storage of accept-reject matrix 
Counter - no. of sieves rejecting a batch 
Titte "for ‘output 

Trele *for"output 

LLere sor output 

Trove -tor outpur 

Title for output 


Counter - tolerance accepts, rejects and 
percent rejects 


Counter - uniformity and difference accepts, 
rejects and percent rejects 


The State's upper job mix formula band 


Single Variables 


Number of standard deviations adjustment varies 
Value of adjustment of percent retained 

Mean percent retained adjustment 

Number of batches to be simulated (same as NBATCH) 
Number of hot bins (same as NBIN) 


A blank signifies an agcept in accept-reject 
matrix 


Program Symbol 


DIFF 


DIFRET 


DRFN 

DRF'RN 
DRIFT 
PLRo 


HIT@L 


IN 
IPBIN1 
IPBIN2 
LPT 


IPMRT 


IPRINT 


IR 


Meaning of Symbol 
Single Variables 


Difference between retained values 


Difference between sum percent retained and 
100 percent 


Signifies the second drift "2ND" 
Seed eocecandon numbers that determine drift 
Value of drift applied 

Stabting sbacchwofeL,lLnstednl at 

Calculated upper tolerance limit 

SUupSscrapt =. genera! 

Counter -— number of drifts 

SUDSCHIDE = bins 

Temporary value of the primary sieve number 
Counter indicating if plots are required 
Counter indicating if plots are required 
Counter — drift 


Number of tolerance rejects 


Counter — correction on uniformity difference 


Counter — drift 

Subscript =— general 

Input reader = 5 

Pirst bin<with primary sieve 
Last bin with primary sieve 
Current plant number 


Rounded faxed value of the sum of plant.mix 
ratios 


Number of batches to be printed 


Seed for next generated random number 


Program Symbol Meaning of Symbol 
Single Variables 
ISAM Counter - samples 
LISAMP Subscript - sample 
ISIEVE Counter — sieves 


ISTARI Temporary Gritt istart value, first drift 


ISTAR2 temporary; drift start value, second drift 
ITEST Switch value 

ti bod) Counter -— steps 

Tool Counter — steps 

J Subscript - general 

JMF' Job mix formula 

JSIEVE Sieve number of primary size 

K Subscript - general 

L ~uUDSCript — general 

LOWT OL Caiculated “lower tolerance’ limit 

M Subscript = general 

MARK Counter indicating af JMF is. outside State band 
MBIN Subscript - bins 

MINC Difference between samples 

MSIEVE Counter -— sieves 

N subscript — genera | 

NAME Name of a particular sieve 

NBATCH See BATCH 

NBIN Subscript - bins (also see BIN) 

NDRIF Sweten andicatingvifedriftois desired 


N@PLT Number of plants to he simulated 


Program Symbol Meaning of Symbol 


Single Variables 


N@PRIM Number of bins having a primary sieve 

N@®@RRN Seed for normal random number generation 

NP Subscript, — <disStrtoutllon, come 

NPT Number of points: in a adLisSerioution 

NSAMP Total number of batches simulated 

NSIEVE Number of sieves (same as SAMP) 

NSTEP Switch indicating different sampling interval 


(same as’ STEP) 


NTSAMP Total number of samples 

NUMAC1 Counter - uniformity accepts 

NUMAC2 Counter - uniformity difference accepts 

OUT Output printer = 3 

PASS2 Temporary value of percent passing 

PCREL Mean percent retained value 

PCRRN Seed for random numbers that determine RETAIN 
values 

PMRI Plant mix ratio values times 100 

PMRT Total sot vol anima tL 

PRNT Number of batches to be printed 

* R Signifies a reject "REJ”" in accept-reject matrix 

RET Temporary value of percent retained 

RETAJ Adjustment of retained to total 190 percent 

RND A random number 

RNPCR Percent retained random number 

SAMP See NSAMP 


SAMP2 Temporary value of the total number of samples. 


Program Symbol 


SDIF 
SIEVE 


SIGMA 


SPRIM 
STATEH 
STATEL 
STDPAS 
STDST 
sTDT 
STEP 
SUMPAS 
SUMPSQ 
SUMRET 
aol: 
homie 
TESTDF 
TESTUN 
VALADJ 


VALSTD 


SUBR@UTINE N@RMAL 


Meaning of Symbol 


Single Variables 


Imposed uniformity difference 


VE ee 


see NSIEVE 


Number 
varies 


of standard deviations percent retain 


Imposed uniformity limit 


Imposed upper tolerance limit 


Imposed lower tolerance limit 


Standard deviation of percent passings 


Imposed two sigma tolerance limit 


Two sigma limit of percent passing 


see) NSTEP 


Sum of 


Sum of 


Sum of 


Number 


Number 


Number 


Number 


percent passing 

Squared percent passings 

the percent retained values 

of sieves rejecting -— tolerance 

of samples rejected by uniformity difference 
of bins rejecting - uniformity difference 


of bins rejecting - uniformity 


Mean standard deviation of adjustment 


Mean standard deviation of 


percent retained 


Arrays 
Cumulative normal frequency values 
Number of standard deviations 


Program Symbol Meaning of Symbol 
Single Variables 
FX Random number generated (from main program) 
I Subscript - general 
I@ Number of the output printer 
J subscript i— general 
NPT Number of points in cumulative normal distribution 
+ VNQ@R Number of standard deviations 


SUBROUTINE GETVAL 


Arrays 
FN Equivalent of ALLRET main program 
W Bracketing mean adjustments of percent retains 
WwW Bracketing standard deviations of mean adjustments 
X Bracketing random numbers 
¥ Bracketing mean percent retains 
Z Bracketing standard deviations of percent retains 


Single Variables 


+ AMEAN Mean percent retained 
+ AJMEAN Mean adjustment of percent retained 
FACTOR Interpolation factor 
E Subscr ppt dsigenerak 
I@ Number of output printer 
J subscript — general 
K Subscript - XN value in FN 
KK Subscript -—- AMEAN value in FN 


KKK Subscript — SIGMA value in FN 


Program Symbol 


K4 


SUBRGYUTINE GRAPH 


Meaning of Symbol 


Single Variables 
Subscript - AJMEAN value in FN 


Ssubscript..— SIGAJ value inv FN 

Sieve, number 

Bin number 

Number of points pm distribution 
etandard deviation adjustment 
standard deviation percent retained 


Random number 


Arrays 

Number of rejects 
Preguency storage 
Dashes tor output 
Counter for plotting 
Tie le, Ono Loc 

mingle Varleables 

An "#" 

Ag bud ak 

Maximum frequency 
subse riota — genera. 
Staring factor for splet 
Scale Of print out 
USS as general 
SUbSC rips .— general 


Maximum scaled value 


Program Symbol Meaning of Symbol 


Single Variables 


- N Number of sieves or bins 
NUM Number printing spaces available 
QUT Number of the output printer 
SCAL Scaling factor for plots 
X Temporary value 


SUBR@UTINE RANDU 


Single Variables 


- IX Seed for random number 
m4 Seed for next random number 
+ YFL Generated random number 


SUBRGUTINE ZERQUT 
Arrays 


- FREQ Frequency of rejects 


Single Variables 


- wN Number of sieves or bins 





* These variables are specified ina data statement 
+ Values sent back to the main program from subroutine 


- Values sent to the subroutine from the main program 


APPENDIX D 


INPUT CARD SPECIFICATIONS 


Card Column Variable Field 
No. Numbers Description of Variables Name Specification 


Card Number One 


it: Number of Plants NZPLT BL 


Card Number Two 


1-5 Percent retained random no. seed IRN (1) 15 
6-10 Normal random number seed IRN (2) 5 
11-15 Drift random number seed IRN (3) I5 


Card Number Three 


5 Tolerance 1/2" sieve STATE (1) P52 
6-1) S 1/4" sieve SPATE (2) F5.0 
11-15 i 1/8" sieve STATE (3) F5.0 
16-20 s #20 sieve STATE (4) Bo 0 
at-25 4 #80 sieve STATE (5) F5.0 
26-30 B #200 sieve STATE (6) 5.0 
Card Number Four 
Che Upper JMF Limit 1/2" UPLIM(1) F5.0 
6-10 ed : tye UPLIM (2) PS. 0 
11-15 ‘ 3 i 1/48" UPLIM (3) F5.0 
16-20 r " ’ #20 UPLIM (4) F5 0 
21-25 : #80 UPLIM (5) . F5.0 
26-30 a « . #200 UPLIM (6) PS 
Card Number Five 
I=5 Lower JMF limit 1/2" LOLIM (1) F5.0 
6-10 W uw W ae L@LIM (2) F5.0 
11-15 . C ng » by ono ity LYLIM (3) age sf. 
16-20 Es a x #20 LOZLIM (4) FS. 0 
21-25 : " x #80 L@LIM (5) F5.0 


26-30 * 2 #200 LOLIM (6) FS 


p=2 


Card Column Variable Field 
No. Numbers Description of Variables Name Speci fication 


Card Number Six 


a> “Un@Pornicy LoeLe SPRIM F5.0 
6-10 Uniformity difference limit SDIF P5770 


Card Number Seven 


1-5 Number of batchés ; BATCH FSG 
6-10 Number of bins BIN F5o0 
11-15 Number of sieves SIEVE F5.0 
L6e-20 “Drettaactuator DRIFT P20 
21-25 Number of steps STEP 1s rf 
26-30 Number of batches printed PRNT Fon 


Card Number Seven A (Optional-insert card only if STEP=0 on Card Seven) 


1=5 ValuéfotyStep STEPS (1) F5<0 
6-75 Up to 14 more step values TOL STEPS iS!) 14F5.0 


Card Number Eight 


i=-5 a PMR ieee. PMR (1) 5.0 
6-10 PMR Bin #2 PMR (2) le ape) 
L1=15 (PMRABIRATS PMR (3) Bo <0 


Card Number Nine 


1-5 Primary Sieve Bin #1 PRMBIN (1) F5.0 
6-10 ks "s Bin #2 PRMBIN (2) B5.0 
tril ES . o Bin #3 PRMBIN (3) E530 


Card Number Ten and Eleven, etc. 


A pair of input cards are used to describe each sieve of each 
bin. For example, a three bin, six sieve plant requires 18 differ- 
ent two-card sets. 


The first card of each set defines the bin number and sieve 
number. The second card describes the distributions of the mean 
percent retained, the standard deviation of the percent retained, 
the mean adjustment, and the standard deviation of the adjustment 
for up to three points on the distribution. A distribution may be 
described by up to 40 points (39 assumed straight line segments by 
additional cards placed behind the second card. 


Figure D-l shows the distribution curve for the percent re- 
tained versus the cumulative frequency for the 1/4" sieve of Bin 
#1. The curve in this example is approximated by two line seg- 
ments and is defined by three points. The following page shows 
the procedure for placing this data on a two-card set. 
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Card Column Description Example Variable Field 
No. Numbers of Variable Value Name specifi cateon 


Farst (Gard of.ceat 


Bin number 15 IBIN 1a 
2-3 Sieve number 2 ISIEVE 12 
4-5 Number of points 3 NPT Le 
6-9 Sieve name 1/4" NAME A4 


Second Card of Set 


1-5 Cum. frequency Oe ALLRET(1, 2,1) Pow 
6-10 Mean % retained 65%. ALLRET(1,/2, 2) ES.0 
LInUS We Stdeudev.eeorets 1.0) SAYRE T (en 203” 5.0 
16-20 Mean adjustment Li2s 7 ALLRET e254) BS 8G 
21-25)  Std.devdvead;. 14 | RUDRerd, 2.5) F5.0 
26-30 Cum. frequency 0.2) (-ARURE fie, 2,6) Boe 
31-35 Mean % retained 4s ALLRED 1,2 ,:/) Ebi 
36-40 Std. dev. % ret. lees ALLRET(1, 2,8) Bo. 
41-45 Mean adjustment 0. 6°") [ABD RE La, 2504 F0 G 8, 
A46-50| Std. dev. adjé 0.8 .ALDRET CL, 2 S10} F5.0 
51-55 Cum. frequency Ve0  PALGRET ie 2 oS eae) 
56-60 Mean % retained Sl. ALLRETCL, 254) 2) OO 
61-65 Std. dev. % ret. 2-0 ADURE GI yhe yy Bee 
66-70 Mean adjustment -1.1 ALLRET(1, 2,14) Fa. 


71-75) Stdf# dev. adj. 1.3) WALLREDGI ee) F5.0 
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APPENDIX F 


LISTING OF FORTRAN PROGRAM 


ASPHALT PLANT SIMULATICN 


MAXIPUM NUMEER © F BATCEES NOW AVAILABLE IS 500. 


HERR EKKH HANNAH RAKSHA AMAA ADHD ARH HKHARKRHHHHHH RHE RRR RES 
PERINITICGA Ch ARRAYS AND VECTERS 


AU) FIA RRAY BYE IN ShYe SPBVE) BY TOTAL NOG. OF POINTS—A, DISERETE 
CISTRIPUTICN CONTAINING AS POINTS THE PERCENT RETAINED 
MEANS A®D STANDARD DEVIATIONS FOR EACH BIN AND EACH SIEVE 
AND THE MEANS AND STANDARD NEVIATIONS FRO THE ADJUSTMENT 
FULF TC ERRORS IN SAMPLING, SPLITTING, TESTING, ETC., TO 
PE MANE TQ THE PERCENT RETAINED. 


PCRET SARR AY RY SRINR EY STEVE=-MEAN .OF THE PERCENT .RETAINED 
FCR THAT RIS FOR THAT STEVE FOR THAT PLANT 
VALSTD -—-ARRAY CF STANCARC DEVIATICNS CORRESPONING TO PCRET 


ACJMEAN -ARRAY DY BIN BY STEVE--MEAN OF THE TESTING ERROR 
DLE TC SAMPLING,SPLITTING,sTFSTING,ETC. 


VALALJ —-ARPAY CF STANDARD DEVIATIONS CORRESPONDING TO ADJMEAN 


NPTS -ARRAY BY RAIN RY SITEVE-NO. CF POINTS OF THE ALLRET OR 
ACJRET OISTRIBUTICNS. 

SNARE —-VECTCR PY STEVE-THE STEVE NAMES IN A4 FORMAT. 

PPR -VECTCR BY BIN-THE PLANT MIX RATIOS IN INCREASING BIN 
SEQUENCE. 


PRMBIN -VECTCR BY BIN-THE VALUE IS EQUAL TO THE PRIMARY SIZE 
SIZE STEVE FOR THE BIN WAS SPECIFIED. 


RETAIN -—ARRAY BY BATCH BY BIN BY STEVE-THE PERCENT RETAINED FOR 
THAT BATCH FOR THAT BIN FOR THAT SIEVE. 


PASS -ARRAY BY BATCH BY STEVE-THE PERCENT PASSING FOR THAT 
BATCH BY STEVE. 

STATE —~VECTCR BY STEVE-THE STATE*S IMPOSED TOLERANCE. 

UPLIF -VECTCR BY STEVE-THE STATE'S IMPOSED UPPER LIMIT. 

LOLI -VECTCR BY STEVE-THE STATE'S IMPOSED LOWER LIMITS. 


MSTEPS -VECTCR RY NO. OF SAMPLING PROCEDURES DESIRED-EACH MSTEP 
CONTAINS THE INCREMENT ON WHICH TO SAMPLE. 


ty 
L] 
nN 


XGY CO Sere eG Canc) 


Ga 


PRU RBWNHRK AAA 


CEShet -~VECTCR 2Y NC. CF RATCPES-IF THE VALUE IS 7E8Q, THE PATCH 
wAS ACCEPTED CN TOLERANCE TEST PROCECURES. LF CHE VWAEVE 
(S CTHER TRANGZERGS TT TS. FOUAL) TOs THE MG OR ey 
BATCH WAS REJECTED Ch. 


RHR MRE RRR RRR K RAAB RHA ERR ARKH RRR AH HERR HERRERA RHEE 


CCUALE PRECISION PINNAM 
REAL JMF,LOKTCL,LOLIM 


INTEGER CUT, PCRRNySNAME( IG) sTESTIe TEST2, TEFTUN,TESTOF,TESTOL, 
LERF RA pCUR(Z) 


\ 
PIMENSTON ALLRET( 200) ¢IRN(4) ,NOTS(S,10),1START(2),y 
IRNO?) pBINNAM(5) ZRETAIN(SECO, 5,19) sPMRID) 
2PASS (500,10) »PRYBIN(5) sTESTCL(5C9) TEMP (2) pPARKY( 500,10), 
ZIFLACC (3) »URACC(293),BACC(2,2) ,ACGREJI 2, 2),REJACC( 252) yERE I 25 7)% 
ASTEPS(15) pFREGEV) gFREGL(19) ,FREQC2(1C6),FREQ3(10),STATE(1C), 
STV (2a STITT 220) TUT S(20) TET 4 (207571150207 UePCIP Clay, Celie tion, 
EMSTEPS(15) pPREMAX (2) pLREN(2) pPORET (5919) pVALSTON(5,19) yASMEAN( Sy 1 
TEV oVALANI(S,19) pg JME (19) sSTATELUIC) SIATER( IU) 


CCNWCK OUTZALCLRET 
EQUTVALFEACE (CERNCL) sNCRRN) pC ERN(C2) »PCRRN) » CIRN( 3) sDRFRN) 


PATA RENNAM/*DIN 1L*, "SIN 2¢,°DIN'3%," [MH 4%, "GIN S°/, 
LPLANK aR oP RENS 9 fg t*RESt,F1ST#,"2NCI/ . 
CATA TITL/*NG. "e'CF Stat TEVE'y*S A *,*SAMP',°LE [*,°S REt, "JECT, 
LIED CHEN 81a FPS, TE L2/INT. 8 PORT Set IMES gts EAC Ryo ttetias 
2'EVE *,*REJE*,"CTEDS,® SAM*,"PLES', 108" "/,TIT3/'°NO. 'y'OF Bty 

ZENS 49A SAtSIMPLE',* TS t,*RFJE%s *CTED', '-UNE', TFORM*Z "ITY %y 
AQa® ©/, TIT4/°TCTA,*L NG ,*. FEY%,! REIT, "ECTS", © BY *,* EACH", 
GUERIN DEZEO8 FS TITOS NGS 8, "OF RUSTING (OY AOSAM SMpien, | 1S as 
EMREJE*SOCTED sg *-UNT "gt FORMS tITY "gD FFT, (ERENT ICE thet ', 


FCO2MsaTs #CR £70 LISTS FOLLCW 


FCORMAT(415) 
FERMAT(5F5.0) 
FORMAT(T1L,212,44) 
FORMAT (14F5.C) 
FCRMAT(T2) 
FCRMAT(T1) 
FORMAT(///1CX_,*PLANT MIX RATIC FOR BIN NOQW*,12,' IS ZERO OR NEGATI 
1VE.*/) 

FORMAT(CLCX,*S EMULATION OF PLANT NOQ.*,12,° WILL TERMINATE DUF TO ER 
1RCR IN INPUT.) 

FCRMAT(///1CX,*THE SUM OF THE PLANT MIX RATICS IS NOT EQUAL TO 100 
1 PERCENT 227) 

FORMAT(S//20K%y*PLANT MIX RATIOS*/20X%_,"BIN',5X, "PERCENT! /) 

FCRMAT(19X,13,7X,F5.21) 

FORMAT(//26X%y"PERCENT PASSING*,38&8X,"PERCENT RETAINED" //37X_ 
L*PERCENT*/10X,* BATCH" py L457X_*STEVE® y6X_ "PASSING 1 7X, "SIEVE', 
25(9K,A5)) 

FORMAT (27 XgA4 gS XgFle2y LIX pA4yS(TXyFTe2)) 


FCRMAICLHEL//42X,'A S PHALT PLANT SIMULATILION*S// 
11CX,"*NUMPER CF PLANTS TO BE SIMULATED IS*,I2,'.'*) 
FCORMAT(///20X% 940 "*NO] OF *45X)/20X9 *BATCHES' y5X_y "BINS" y6Xy*SIEVES', 
LOEXe*STEPS*/21X e914 97TXy3(13e98X)/S/) 
FORMAT(LCX AS 99X yp Fl e2e6( 8X eFTe2)) 

FORMAT(//39X,*CALCULATED JOB MIX FORMULA AND IMPOSED JMF TOLERANCE 
1S  *.//23X%y4("*CALCULATED® 95X) 9 1X9 30" IMPOSED" y 8X)/10X_ "STEVE" ,11Xy! 
2IMF%94X%—2(5X_*TCLERANCE! 54X_,"LCWER LIMIT',4X,"UPPER LIMIT")/) 


FOE MATOUS/2X "PRINT RESULTS CM CETERMINING TF FINAL PERCENT PASSING 
1 fF A SAMPLE FCP CACK STEVE FALLS WITHIN TOLFRANCE LIMITS FOR THE 
PSTEVESPSIOX,*LF=REJ, A REJECT HAS PEEN OBSERVED CP //5X_"SAMPLEN/S TX, 
SINE o Pipl 2A STEVE! 14 

FERMAT(SX,1S 99X91 (A3,9X)) 

FORMAT(//ICX, "TEST GN UNTFGR*ITY GF SAMPLES WITH PRIMARY STZCL PEPE 
LENT. PASSANGE" FO. 3/S 10X52 SAMPLT NOs. 1).502%,85)) 

FFAMATILCX I 5 eIX%e5( 834K) ) 

FOIMAT(S/LCX eg STEST CM DLIFFEQSENCE BETWEEN SAMPLES*/1LO%, "DT IFELRENCE 
IMUST BE ECUAL CR LESS THAN ,FG63—%e"/S/1LOX— "SAMPLE NOS 04 952%, 
245)) 

FCORMAT(9OX, 14,9 '-',1496%,5(43,4XK)) 

FORMAT(//1CX%_"*KC PRIMARY STZ SPECIFIEN FOR. ANY BING //10X, 
L°UNTFORMITY ANC TCLERANCE CLEFERFNCE ARE NOT ATTEMPTED. YS) 

FOIMAT(/2CX,* RANDOM NUMPER STEEDS USFT FCP THIS PLANT TO INITIALIZE 
L RANTCH NUMRER CENERATORS ARE—"/20%,"NORRN=",T1Cy"%, PCRRN=",110, 
28 en PRERNE?' TC) 

FORMAT(IE1L,1GX,"IME FOR THE, A4,*STFEVE TS GUTSICE STATE LIMITS®) 

FORMAT{(LCX,"THE SIMULATION OF PLANT NO2*%,12,* IS TERMINATED SECAUS 
1F THE JMF IS CUTSTINE THE STATE LIMITS.*//) 


2aG EReve a is fade er MP EE - OF PF E S Ue FP Sha) 


3201 


TORMNAT(L2X,52(¢ #*)) 
FORMAT(LGX sf ee, L LX, 3098, 7K) 2X e9( HY, TK) ) 

FORMAT (1LH+ 532% ,4(°NC. CE°Z19X)) 

FORMAT(VH+,26X,"ACCEPTS REJECTS REJECTED", 20" ACC-ACC*'),2(" REJAR 
LES") 920" ACC-BREIN) -20% REJ-ACC")) 

FORMAT (VE +,15X_,"TCLERANCE® ) 

FORMAT (LE+,14X%,"UNTFORMITY®) 

FOPMAT(1H4,14X,"CLIFFERENCE® ) 

FORVATULEF p22 XFS og 3X gf SoD PX gp FT ola lX 9 4( 2X pF 50g 2X pF Te2)) 

FORMAT (LB4 pp LEX, "TEST 9 SK%e20"NDO. CEP2X), "PERCENT EVENTS PERCENT? 
193(" EVENTS PERCENT!®)) 

FORMAT(IE1L,LOX,* THE FREQUENCY PLOT ENTITLED *,2UA4/17X%_"WAS NOT 

1 PLETTEN FCR ALL FREQUENCY CLASSES WORE Z7FEROG.'/) 
FORMAT(2CX,*SAMPLES ARE CRETAINED FROM EVERY BATCH.YS/) 
FORMAT(PCKX st INTERVALS AT WEICH A SAMPLE IS ORTAINED-'/ 

PIS C27 PPI7VS*, ACH VEVERY. TH PATCH.'/)) 

FORMAT (//51¥%¥ 9 *IMPCSED*"/20X,* STEVE" s4X%,"*STATE IMF LIMITS' 4X, 

L' TOLERANCE */1C(21LX eA eg SX gF 5009! —" pF 5 eC yp BX, FGel/) ) 

FORMAT {2CX,"THE UNIFORMITY LIMIT IS",F6.2,"PERCENT AND THE UNIFOR 
IMITY DYFFERENCE [S*,F6.2,*PERCENT® ) 

FORMAT (23XyA5_92X_A4) 

FCRMAT(/20X%,"°AC ERIFT SPECIFIEC FOR THIS PLANT.") 
FCRMAT(/2CX,"*CRIFT SPECIFIED FCR THIS PLANT.!') 

FORMAT(3CX,*START CF EURATICN MAX.CRIFT'/31X,"ORIFT IN BATCHES 
1 (PEPCENT)*) 

FORMAT(2CX,A3," CRIFT el S plies l Sy. Me F Ga 1) 

FCORMAT(2CX,*CRIFT wILl BF AMPFC TO TRE PRIMARY SIZE STEVES ON THE 
LFCLLOWING BINS-—*/25X,"BIN 1°) 

FORMAT (25X%,45) 

FCORMAT(2CX,*AKCD THE SAME AMOUNT OF PRIFT WILL BE SUBTRACTED FROM T 
LHE PRIMARY SIZE STEVE ON BIN 2.°) 

FORMAT (2CX_e"THE PRIMARY SIZE STEVES RY BIN ARE—'/29X%,"STEVE") 
FORMAT(/1GXy,*SAMPLELD EVERY',14,°TH BATCH.'/) 
FORMAT(/1ICX,*ALL* 9149" BATCHES FXAMINEN.'/) 

FORMAT(/10X,"ANALYSIS OF EACH CF THF',14,* BATCHES FOR THIS PLANT 
17S AS FOLLCKS-"//) 


1CCCC FCRMAT(LEL/S/////) 
ECC GL FERMAT(4G5SX,*#* 40Xy te! ) 
1CCC2 FORMAT(45X_,"#ASPHALT PLANT SIMULATION EXITED NORMALLY*#*) 


1CCC3 FORMAT(45X%942(%#")) 


LLCCC FCRMAT(1H1,30%,'A S PH ALT PL 
"4! 


C 
C 


AWN, Ta Seek My eke A Tok OA Oa 
) 


Datu bekEoN UF DUFF ION, * PLANT NOL." F2 


ie AY sim JOS op Wes a J ead ge Rs YS sae Sai en 


oO 


CEFINE T/C UNITS 
[N=5 
cunts 


ACPLT-NUMPER CF PLANTS TC PE STMULATED. ITF, GNLY 1 PLT. ASSUMED 
READ (IN.5) ANCPLT 


CEFING -ROKECR NUMBERS 

IRA-RANDCM NUMBER SEEDS. (MAX. CF 4 AVAILABLE) 
REAC IN RAANFCM NUMFER SEFCS 

RFAR (CEN, LI IRN 

ACRRA=TRN (1) ; 

PCRRA=IRA(2) 

CREN=TRN(3) 

REAC(IN,4) STATE 

READ(CIN,4)UPLIM 

REAC(IN,4)LCLIM 

REAN- UNTFGRFITY AND UNTFERMITY SLPFERENCE LIMITS 
REANCIN,2)SPRIM,SCIF 

TE “CREPE T PL 2 eZ aes 

ACPLT=1 

PRINT HEALING 

WRITF(CUT,LCS)NCPLT 

RE VYLDL FRCP _eXCPET 

PRINT KEARIANG FCR PLANT 

WRITE(CUT,LLOCU)IPLT 


READ PARAMETERS 


PATCH=TOTAL NC. CF STIMULATED BATCHES REQUESTED 

BIN=NO. CF HCT BINS 

STEVF=N0. CF STEVES 

CRIFI-IF Cy, NC DRIFVICESTRED, OTHERWISE OR IRM OINGLUD ERIN 
DETERMEINATICN CF PERCENT RETAINED. 

STEP=NC. OCF STEPS, “FPeEs INEREMENT GN “BATCHOSOST RATVEDR “EVERY 
SIMULATED BATCH AS PENIFIEC RY “STEP FA SIMULA TEDYSAMPLE . PS ""RTAENED. 
PRAT=C ITF ANALYSIS CN ALL BATCHFS IS TO Bt PRINTED. 


READING 4)PATCHyPINgSTEVE,CRIFT,STEPgPRNT 
INITYALETZE TCTAL NC. PATCHES 

NEATCH=BATCEH 

NTSAMP=RATCH 


NPTN=BIN 

NSTEVE=SITEVE 

NCRIFH=ARS(CRIFT) 

NSTEP=STEP 

IPRNT=PRAT 
WRITE(CUT,210)NCRRN»gPCRRNyCRERN 
WRITE(OUT,LIGINBATCHsNCINsNSTEVESNSTEP 
TF(NSTEP)14,14,12 


SET NSTEP=1 IF ZERO CR NFGATIVE 


READCIN,4)(STEPS(I),IT=lyNSTEP) 

Cf 6666 T=lyNSTEP 

MSTEPS(T)=STEPS(I) 

READ PLANT MIX RATICS BY INCREASING RIN NOS. 
READ(INs2) (PMR(J) sJ=1yNBIN) 

REAN WHICH BINS HAVE THE PRIMARY STZES. VALUE IS EQUAL TO 
WHICH STEVE IN DESCENDING SEQUENCE IT IS. 

READ(IN, 2) (PRMBIN(MBIN) »MBIN=1,NBIN) 

TIPRINL=0 

IPBIN2=0 
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i oe te ae | 


WEE Cty OV TA 


ie see ce Ne NSD he Lis 

TF CPRYBIAN( TP IN) 222292222974 
FIAT PFEGIANERG ANE ENDENG CINS WHICH HAVE PRIMARY 
S177 STEVE "EC hehe 

IF (FPOPENTI24,34,35 

TPE TSP =TeE TA 

VER INZ=T EEN 

CENT WHE 

TR (NSTEP)S3,52494 
WVRETE(CUT,4¢C1) 

CH TS, s 

WRITE (CUT ,4CO7) (MSTEPS(T) » F=lLeNSTEP) 
WRITE CCUT,1C3) 

CPR Ph res 

PYRY=o 4 

TA 822 T= Le Work 

PRE T=PKR(TPeLCC.S 
WRITE(RCUT, LOS) Lore ot 

LF CPV ACE Vs eb 34.22 
WRETECCUT,LCO)T 
WVREFECCUT,LCLILPLT 

er OU. SA Te 

PYRT=PKRTEPMR [ 

TPVRT=PMRT+C.5 

[IF CTPMATLEC. LOC) GC Te 15 
WRITC(LUT,1IC2) 
RPITECCUP Ss LELIIPUT 

GC pO SET 


PEAT ALLRET PISTRIRUTECNSC ANY CRHER AS LONG AS REAN AS AMSET hs 


EC) 32) T=), hOEN 

GR 32 J=1 sNSTEVE 

READ CINe SPIF ING ISTEVE Ss NPT yNAME 
GET cSILVE -SEZE “NAEYE 
SUAMECTSEEVELENAME 

NPT=A9T#5 


NPTS(TRIN,ISTEVEVSNPT 

REAM (CIN,4) CALLRET OND) SNP=L NPT) 
CALL RANEU(POCRRNg TReRNPCR) 
PCARA=IR 


CETATN MEAK ANE STANRARC CEVIATICN FOR PERCENT RETAINED AND 
MEAN ANAND STANCARR CEVIATICN FGR ADJUSTMENT FACTOR DUES TO ERRORS 
IN SAMPLING, SPLITTING, ANE TESTING. 


CALL GETVAL (NPT »RNPCRePCRET( TP ING ESTEVE) pVALSTCCIBIN, ISTEVE) sAIME 
LAN (TRIN, ESTEVE) pVALACICIE INS STEVE) ) 

WRITE (CUT 493) (SNAMECTSTEVE) sLOLIM(ISTEVE) SUPLIMCISTEVE) » STATE 
LUISTEVE) p ISTEVE=1yNSTEVE) 

CO 8f88 [=1,2 

CUR(T)=9 

[F (TPBUNL) LLG, 1169115 

WRITE (NUT ,404) SPRIM,SCIF 

WRITE (9UT,420) 

CO $599 [=IPEIN1L,IPBIN2 

IC=PRWBIN(T) 

WRETE(CUT, 405) GC INNAM(T) » SNAMEGIC) 

CRIFT SECTICA(SETUP)-CNLY APPLIFS FOR-1A TOP WITH 3 HOT BINS. 
ITHRU=C : 

IF (NCRIFI116,1165117 

WRITE(CUT,406) - 

NORIF=0 

TSTARI=NBATCH4I 

TAMCN=0 
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7777 


134 


(Cue Cai 

[TF (ITTHRL)54,123,54 

NOR TE=1 

CALL RANRUE(RRERN,TR,RNC) 

IE CRNO=C 75) 1a hes 19 

NOR TF=2 

JANCN=ANDRTE 

ISTAR2=0 

WRETE COUT 497) 

WRETE(CUT,4CR) 

Rie MANE RATE, ADA SAE Re Ae 

CALL RANOL( IR CRERNeEIRST) 

BIR SIAC. URES Vee Fe 
ISTARTCICFT)=HFIRST#+LISTAR2 

CALL RANCU(CRERN, TR a RNC) 

CALL RANCUCIR LRFERA,OCREMAX(ICFT) ) 
RBNG=6 So SRAE 85C CG 

CUR CTINET )=254+kRNG 
PREVAX(LLET I H3.C44.O8#NREMAX(ISFT) 
ISTARZ=ISTARTCIOETI40UR(CICEFT) 
[R=aDRFRN 

WRITECCUT, 499 CRENCICET) g ISTARTCICET) »DUR( LOFT) »sONRFMAX( IDET) 
PREMAX(TCFT)=9.8#fRERAX( IDET) /(2.08CURCIGFT)) 
TSTARL=TSTART(1) 

CALL RANCULERERN, TR gRNC) 

NCRIF=1 

WRITE(CUT,410) 

IF (RND—6.56)121,120,122 
WRTTA(CUT,412) 

NCRIF==2 

DRFRA=ER 

Geis ple 4 

CALL RANCUCTR BDRERNeaPND) 

TE CRND=H=C S01 123, 1230124 
NORTFANNRTF AL 

TF (RND-0.89)126,126,125 
ANCRIFHANDRITF +1 

WRETEC( CUT, 411) °R INNAM(1), f=2,N9R TF) 
ICFT=3 

Nh Ey ples 

WRITECOUT,LIOCO)TPLT 

TF (TPRNT) 149,149,150 
WRITE(CUT,SC2)NRATCH 


THIS SECTICK CBTAINS PERCENT RETAINED CFOR GACH “ShEVE, 
INCREMENTS CN STEVES», PINS» AND SAMPLES. 


TC 66 [SAMP=1,NRATCH 
IF (TTHRL) 136,129,136 


CHECK FOR CRIFT CN OR NOT. 


IF (LSAMP-ISTAR1)129,128,128 
IF (LAMON-ICFT)127,1279139 
IGC=1G0+#1 

ISTAR2=1STAPI-1 

GF TC (131,1329133)s,1GC 
ICFT=1 

GC TC 134 

TSTARL=TSTART(2) 

INFT=0 

GC FC 136 

iCFT=2 
TSTAPL=ISTARL4#0UR(ICFT) 

Gc TC 129 


oO 


it 


) 


NOP RT =F 

[er f=5 

PETAL Ss PATCH 4lL 

Cr o77 FPiked RATS 

PICKLP PQTKMARY STZt STEVE Fe THIS PIN TF THER IS UNE. 
JSTEVESPRMOIN( TRIN) 


SET SUP CF PERCENTESGE PE RAIMETR VALUFS FO TERMS) (BALLS FS USED WATER 


1 CY RESET PERCENTAGE RETAINED FOR & RIN TO Loo PERCENT. 


S[veryen er 
fe Pe PSPUVE =P SU VE 


PICK bP 2ONFCE KUEMEERS 


CALL RANETLOINEPRA,ST8® sRN(1)) 

CALL @ANPUCTRgNCERNSENGZ)) 

CALL PERMSL ERA ( 1) oSTGMA) 

ee Oe AAT Se EME PAPER” 

FALL NSEPMALERAN(2),ACISIG) 

RE TATE (TSAME SEP IN ESTEVE) SPCRET(ICIN,ISLEYS J+SICMASVALSTICUISIN, 
ENE } 


Gt FSIS PP) oe ESET CR PVE EG .PRROR SLOTS! SAMBLAINGs,) SPLATTUING AND 
TEeT PRE. 
SCJUSTINIMECN (TRIMS LSTEVE DFAT ISTIGRVALSTO(ISIN, ISIEVE) 


GET PEICERTAG® RCTAINEL FER A SAMPLE, 
RETAIN (TSAMFE, TE IR ,ISTEVE D=SFTAIN( TSAMP,ISIN, ISTEVE )4ADJUST 


ACh ceeck TO SEF TF CRLIFT TS TG ®F APPLIED TO THE PRIMARY SIZE 
STEVE C2 ROT, 


Te OFT DS ld PSS 

ee tITFT) 145, 1350185 

PF (SST VE-PSTEVED1L355137,135 

PS CeTEy TC FE ASPPLIEN Th FTES GN. 

TF CTR EIN-TAPSONEE FEY )135 9.136 49.1 35 

PATCH HTS AV F-TSTAR2 

TF SC» MAKE ORIFT CORRECTION CN PERCENTAGE RFETAINEC. 


DREFTHCREMAX CTCF TD) ePATCHe(L.9-(RATCH/NURCIDFT)II/(0.2) 

TF (NORTE 139,135,146 

CRE T=-OR TF T, 

RETAIN(P SAWP, TEIN, ISTEVE )ERETAIN(OISAMP, IRIN, ISTEVE) +ORIFT 


AFTER ALL THIS, CHECK PERCENTAGF RETAINED FOR ZERU OR NEG. 


TF (RETAIN(YS4MP,FPIN, ISTEVE)I16,16,28 
RETAIN(TESAMP, TRIN, ISTEVE =0.0 


ACCUMULATE SUM FF PERCENTAGE RETAINER FOR ONE STEVE SIZE. 
SUMPRET=SUMRETHRETAINOISAMP, IP IN, ISTEVE) 
EXC CF STEVE INCCEMENT. CONTINUE ON BINS. 


RYFRET=1LCC.C-SUMRET 
LO 77 MSTEVE=LANSTEVE 


CCRRECT PFRCENTAGE ROETAINEC SC SUM OF PERCENTAGE RETAINED OVER 
A PIA £S 1CC PERCENT. 
RETAJ=H(RETAINGISANP,IBINSMSTEVE)/SUMRET) #DITFRET 
RETAIN(TSAMP,TPIS sMSTEVEJ=RETAIN(ISAMP, [RIN SMSTEVE) 4R Ebay 
CCNTINLE 


IS! 


C ks . 
C EXT £F SIN INCREMENT. CONTINUE CN FORVSATCHES. 
C -o 
-C SET UP FCR FINAL PERCENT PASSING FOR SIEVES. 
C - 
C SET PERCENT PASSIAG IN FIRST STEVE TO (100 PERCENT. 
DASSCTSAMP,1)=1C0C.0 
CC LlL MSTEVE=L,ASIEVE 
PASS2=¢.C 
CAP 22 oe PT N= PN EN 
C SUN PERCENT PASSIAG EVER ALL RENS FOR EACH SIEVE. 


Cee PASS2=PASS2+PMR(VMRIND SRETAIN( ES AMP MARIN, YS TEVES) 
RASS CESAME MS TEVEV=PASSCISAMD eS TEVis ANS 7 
{Fo (MSTEVE-RSTEVE LIT SEE Ls TE 

17 PASSC(TSAMPsMSTEVEF1LI=PASSLISAMPMSTEVE ) 


Pela} CCATINUE 
TF CPPRNT) 11461 14,66 
114 WRITECCUP,ZICS)TSAMP, (RINWAM(T) »TH=1L NP IN) 
CO V4 MS VEN Gale ttS LEVEE 
141 WKREITE(CUT,ICEVSNAMECMSTEVE) sPASS(ISAMP MSTEVE), SNAME(MSTFEVE ) 
2(RETAIN(TSAMP SMP IN SSSTFVE) SMF IN=l, NTN) 
é GUN TENU Ss 


re 


é 
C 
C REGIN SECTION CN TOSTING: WHETHER FENAL PIEREENT PASSING GN 

C A OSTFVE FER EACH SAMPLE (GALES WL RS GWE ANUP Pie te Tit Tt) ie Gulemn tae las 
C FINAL PERCENT (PASS TUNG HIN: MRE CSitBNiEl 

C 


WRITE (CUT,LICCO)IPLT 
WRITE(CUT,2C02) 


PICK UP JMFy STCPASS, ANC TOLERANCES GN HIGH AND LOW SIDE OF JMF 
FOR THE SANPLE. 


JMF FS OVER ALL SAMPLES FCR EACH STEVES SIZE BASED Gee I AAC PERCEAYT 
PASS UNG SECRS EAC RSL EVE. 


Co GonG tact elon sC > 


CALLE ZERCUTINSICVES FREC) 
CALL ZERCUT(NSTEVE,FREOL) 
{C1=1 
VARK=0 
SAMNP=NRATCE 
IF (ADRIF)56557,56 
56 IF (ITTHRU)240957,24C 
aT CE 525 ESTEVE=1L I NSTEVE 
SUMPSO=0.C 
SUMP AS=0.0 
CC 666 TSAMP=1,NPATCH 
SUMPAS=SUMPAStPASS(ISAMP,ISTEVE) 
666 SUMP SC= SLMPSC#IPASS (ISAME, ISTEVFE)#PASS(ISAMP, ISTEVE)) 
JMFCISTEVE)=SUMPAS/SAMP 
STC PAS=SUMPSC-—(SUMPAS ®SUMPAS/SAMP ) 
STPPAS=SCRT(STCPAS/(SAMP-1.9)) 
STC T=2.0#STCPAS 
HITCL=JMF(ISTEVE)+4+STCT 
LCOWICL=JMF(ISTEVE)-—STDT 
STLST=2.C#STATE(ISTIEVE) 
STATEH(ISTEVE)=JMF(ISTIEVE)+STOST 
STATFLC(TSTEVEI=JMFCISTEVE)-STDST 
IF (STATER(ISTEVE) .GT. 100.) STATEK(ISIEVE)=100.6 
IF (STATELCTSIEVE) .LT. 0.0) STATEL(ISTEVE)=0.0 
WRITE(CUT,2C1L)SNAME(ISIEVE) »JMF(ISIEVE),STOT, LOW TOs Riis STATE 
LUISTEVE) » STATEL (ISIEVE) » STATEHUISTEVE) 
TF (JMFCTSTEVFE)-UPLIM(ISLTEVE)) 234,234, 236 
234 TFC JMFCISTIEVE)— rrr ad kel Sy oF 40,240 
236 MARK=MARK4I 


éC 


61 
142 


Ce 555 TEAMP=1,NRATECH 
JSTEVE=ISTEVE 


TFST FOR RICH ANE LOw LI*TTS ON FINAL PERCENT PASSING 
FER FACH SIEVE FCR EACH SAMPLE. 

SESITt WIL CENPATE RO. CF REJECTS: OVER SHEVES 

FER ALE SAMPLES, 


[fF (ISTEVE-1)21921923 
SET Tf TETAL NC. CE PCSSTPLE REJECTS. 


TESTEL CLSAWMP)=ENSTEVE 

IF (PASSCISAMP,ISTEVE)-STATERCTSTEVE))1E,20219 
TF (PASS (CESAPP,ISTEVEV-STATELCISTEVE ) 119,920,266 
PARRY(CESAPMP,TISLEVEI=ER 

RE IFCT 

FREC(ISTEVEV=FRECCISTEVE) 41.0 

IC l=2 

tert o55 

PARRY (TSAMP,TSTEVE)=PRLANK 
TESTCLOESAMEI=TESTCLIITSAMP)-1 

CONTINUE 

WRITE(CUT,LIGSCVEPLT 
wWRITECCUT,263) (SNAME (MSTEVE) oMSTEVE=LSNSTEVE) 
CC 777 TSAMP=HL,NPATCE 

TFATL=TESTCICISAMP) 

TFCQTFAIL) 7779777474 
FRECLTCIFATLI=FRECICIFAILI41.¢ 

WRITE(CUT,204 )ITSAMP, (PARRV(TSAMP, ISTFVE) gs ISTEVE=1 >NSIEVE) 


THIS EXNS SECTICN LTETERMINING TF FINAL PERCENT PASSING FALLS 
WITHIN TTS LCWER AND UPPPER LIMITS. 


PECTS. SECTICR CN TESTING FCR UNTFORMITY AND 
RIFFFRENCE PETWEEN SAMPLES. 


GC TC (75,76),1C1 

CALL GRAPE (NSTEVE,FRECL,TITL) 
CALL GRAFK (NSTFVE,FREC,TIT2) 
CY TF 7 

WRITECOUT,4C0O)TIT1 
WRITE(CUT4CCITIT2 

TF (MARK) 24142741,242 
WPITF(CUT,239)SNAMECJSTEVE) 
wWRITE(CUT,243)0PLT 

IF (TPBIN1)37,37,36 
WRITE(CUT,2C9) 

GE Fe? PERU 

MIAC=1 

NSAMP=NRATCH 

NSTEPS=NSTEPS+41 
NCPRIM=IFPBIN2-I[PBINI+41 

CC 5555 ITESTU=1,NSTEP 
WRITE(CUT, 11000) IPLT 
ITESTI=ITESTU-1 

IF CITESTI)60,61,60 
MIXNC=MSTEPS(ITESTI) 
WRITF(CUT,SOO)MINC 

GC TC 142 
WRITE(CUT,S5SCLINEATCH 
SAMP=NBATCEK/MINC 
WRITE(CUT,2C5)SPRIMs(BINNAM(T) » I=IPBINI, IPBIN2) 
CALL ZERCUT(NCPRIM,FREQ) 

CALL ZERCUT(NCPRIM,FREQL) 
CALL ZERCUT(NCPRIP,FREQ2) 
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GA GU TC i Ge Ps hs beet) 
IC1L=1 
LOg=4 


ERE ANT TEN Gh ARRAY Sa siulcu rieelive 


TERAGE (3) =GUNTAINS TCT4L “TERE RANEG ER WE Ml TNC Gtee TiSiy sit oie 
AKT PERCENT SS JECTED RESPCCTIVEDY, 


THE FOELTRIRG PEETAIN TO BETH UNTECEVITY ANP UNIT PURM ITY Bit Pr eRencr 
VITH TRE “EXCPPTICA TRAT RESULTS? BORSA RERMEEY ORE Ye tick PCA A TINGED 
IN TRE FIRST ROK CP THE ARRAYS “ANB ORESUEPTS = PERT AT ENG 310 
UNUCPGRE LRYe ER PEPERENC E ARE CENTATNEC TN TRE SECON RUW OF THE 


ARRAYS. 
UNACC£ 2,3) =CCRTAINS TETAL YNTEFORMITY ACCEP PShat ee CiiGs 
WieGieN Wo ESV es TNs Ee eles ames 
TC FAL UNTFOERMITY CIFRER EGE AGEEPii are ie otitcolkss 
AKT PERCENT REGECTEC ARE Thee. 
PRCEN Ase) =RESULIS CN EVEMT RETR Ane ACOH Ireis, 
CFURTER TN TST ECLUMN, AND VPEREENTAGE TNZUE Gils 
ACCREI(2s2)) =RESULTS UN EVENT ACCEPT (FIRS Ve lGule suerte tee Nl ual ects 


CCUNTER IN IST COLUMN, AND PERCENTAGE AN aise ia. 
REJACC(I2,s2)  -RESULTS JCNEMEMT REVECT ONS FIRST OCCERT: Tia SEChaG. 

CCUNTER Ik IST COLUPN,- AND PERCENTAGE MING ZA COIL 
RREI(2 927) RIE SUT SIGN EVENT POIs AjeG ois a) © Clete 

COUNTER EN 1ST CPLUMK, AND PERCENTAGE IN 2h COL. 


CTR eT pk ap J (in Yas Fit YB Ye YU a Ja IN Gah oe Sa 0 J ame Ni A ee ge aN aCe a ane) oma) 


TELACT. (3 } se 
VURN © REJECT COUNTERS ANB CTRER GEUNTERS: TQ) ZEees 
Cl: 4444 $ETy 2? 
UNACC(T,3)=C.0 
1CUACE CHIE ert 
TE 44644 J=1,7? 
CRACE( 1,1) =Coe 
PACC(I,J)=C.C 
ACCREJ(I 4,3) =€ .& 
PEUACGIT pweHeoG 
4444 PREJ(IT,J)=C.C 
CC 6£8 TSAMPHL,NSAMPAMINCG 
NUMACI=0 
AUMAC2=0 
TSAN=TSAPP-1 
ZC SS9 TPIN=IPRINLZIPRIN2 
ISTEVE=PRPPINCIPIN) 
RET=RETAIN(ISAMP, IRIN, ISTEVE) 
IF (RET—SPRIWV)25426,26 
2s TEYP(TEIN)=R 
FREC(CIBINI=FREGC(TBIN) 41.0 
Ic l=?2 
Gl Tf 27 
26 TEMP(TRPTN)=FLANK 
NUMACT=NUMAC141 
27 IF (TSAM)28,28,29 
2S CIFF=ABS (RET-RETAIN(ISAM,IPIN, ISTEVE) ) 
TEV (RYIFE-SCIPY Zag 28e2) 
28 NUMAC2=NUMAC241 


PARRY(TSANPyTPIN) =ELANK 
cc Tf 999 

31 PARRY(ISAME,IRINJ=R 
FREGZCTPIN)=FREGZ(IBIN) 41.0 
1C2=? 

SSO. CCATINUE 


a 
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MRETE (CUT ,2CEVTSAMP,(TEMP(MLIN) sMRIN=IPRINI TPE IZ) 
TESTUSTESTICLUCTSAFP) 

TESTUNENCPRIN-NUMACL 

TESTFERNCPRIN-NUMAC? 

CO YC AAC Ss sit ft 

TF (TESTUN) €C,8C,179 

FRECY (TESTUN) =FRECLO(TESTUND 41.2 

Cert, Te 2VP 14 Te? 

Ve" UNO ST CEP? 82 At I 
FRECUTESTEFI=FRECIITESTRE I 41.¢ 

IF (TFST1)2%,28,35 

TCUGCOUT HIT L400 £2341.¢ 

LPPATE REJECT CCUATER FER TOLERANCES 

er aes ees T=]. 2 

MAKE TESTS CA UATFURRITY ANE UNTFORFIVY RIFFERESCE FRR THIS SAMPLE 
SEPT (AC e741 14 1TESI [ 
TESTZ=TESTUN 

Cee 43 

IF (ISAM) @GP oul, 42 

TEST FOR FIRST SAPPLE. SENCE UNTFCQMITY DIFFERENCE HAS MO BEANING 
FOR FIRST SAMPLE, 

TESTZ=TESTICE 

IF (TEST1145,45,46 

MAKE VARICLS TEST AND UPLATE APPRUPRIATE COUNTERS 
SB VACUE VCE FTECO 1S AR SINICA TICON (CE AM ACCEPT... 

TF (TEST29479475428 

tF 1 TEST 2948, 49.,55 

PRCCULIEST, LUSEACCLITESTs 1941.2 

CYENT CF PCTH ACLEPTER. 

ec Tr ene 

EVENT CF REJECT (N TOST2 VARTASLE 

ACCREJ(I TEST, LI=ACCRESCITEST, 1941.0 

oF AT 41 

LPLATE COUNTER CN UNTFERMITY CR UNTFORMITY CLFEFFRENCE. 
EVTNT CF REJECT CN TESTL VARIT4PLE. 
REIACOTITIES Ts LI=REDACCIITEST £1) 41.0 

Ge TC eRe 

FVENT CF PCTHK REJECTEN. 
PREJCITEST,1L)=RREJCITEST,1)I41.0 
UNACCULITEST, 2) =UNACCCITEST, 2941.0 

UPLATE TEST? REJECT COUNTER 

CCOATINUE 

TE (CITESTI)25,91,85 

GG (9G 1094291... 6G) 

CALL GRAPH(ACPRIMsFREQL,TIT3) 

CALL GRAPR(ACPRIM,FREQ,TIT4) 

GOrTC a5 

WRITF(CUT,4°C°)TIT3 

WRITE(CUT,4O0)TIT4 

TFLACC (1) =SAMP-TCLACC(2) 

SAMP2=SAMP 

id Wd Bete Be: Cake ee Bez 


GC TE (52,53) e¢1F1X 

SAWP 2=SAMP2-1.0 

RECREASE AC. ER SAMPLES BY 1 SINCE UNIFORMITY DIFFERENCE HAS ONLY 
A-1 ENTRIES. 
UNACC(LFELX,1)=SAMP2-UNACC(TFIX, 2) 

GET TARLE VALUES 

TCLACC (1)=SAMP-TCLACC(2) 

TOLACC (3 )=10C.C#TCLACC(2)/SAMP 
PPEJ(UIFIX-2)=1CC.C#PRESCIFIX,1)/SAMP2 
REJACC(ITFIX,2)=1CC.CHRESACCIIFIX,1)/SAMP2 
PACC(IFIX,2)=1CC.C#BACC(TFIX,1)/SAMP2 
ACCREJ(CIFIX,2)=LOC.O# ACCREJITFIXs1)/SAMP2 
UNACCCITFIX,3)=1C0.C#UNACC(TFIX,2)/SAMP2 
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2233 CONTINUE 
RKRETE(CUT, 11000) LPLT 
IF (ITFSTL)14341429144 

143. WRETE(CUT,SCLINCATCE 
CC TF 145 

144 = WRETEC(CUT SCO) ENC 

145 - WRITFICUT,207)S0TFs (2INNAM(T)y T=EPCINL» I[PSIN2) 
PC L111 PSAMP=2yNSAMPyMIENC 
TSAM=ESANP-1 

LLLL WRETE( CUT, 208) TSAM, ESAMP, (PARRY(ISAMP, ISIN)» LSIN=IPSINL, EP 21N2) 
IF (TTESTI) 89,924,809 


S2 CRUE Siete Galle f 2 

£6 CAGE GRAPH (NCPRIMSFREQ3S TITS) 
CALL GRAPR(NCPRIM ZFRECZ,TIT4) 
GC Te. 8&9 

&7 WRITE(CUT,4CC) TITS 
WRITECCUT,4C0)TIT4 

B95 WEITTEACUT LLCO PEPLT 


MO 31) 21,4 
Cl TF (312,3135314e315) 91 
312 so TF OC UTESTI)14691469147 
14€ = WRETE(CUT,SCLINBATCH 
CU-TC 142 
147) WRETECCUT, SCC) 
148 wRITE(CUT,3Cc) 
317) =wRITECCUT,3C1) 
ec ¥¢ 311 
313) WRITE(CUT, 303) 
GO TC. 311 
214 = WREITE(CUT,209) 
Ca TOssLy 4" 
215 WRITE (CHT, 3204) 
GC Tc 317 
311 = WRITECNUT, 302) 
CC Sin a3 
J=1-1 
WRITE(CUT,3C2) 
IF (1-1)319,319,32C 
319 WRITE(CUT, 305) 
WRITE(CUT, 308) (TCLACC(K) »K=1493) 
COST es 
32C =o WRITE (CUT, 302) 
WRITE (CUT, 306) 
223 WREFFICUT, 208) (UNACCO eT I) g LI=1)3) 5 (2ACC( SK) Kel 92) 9 (2REI(SyNDy 
IN=197) eC ACCRFE AC SSL) pL=Ly 2) 9 (RESACC( Sy) pM=1y 2) 
IF €1-293225322,321 
321 = WRITE (CUT, 302) 
WRITE(CUL, 307) 
222)  WRETF (CUT, 202) 
318. =wRITECCUT,301) 
§655 CCNTINUF 
IF (NDRIF)58,11111,58 


MINC 


58 IF (TTHRLEILLLI1L» S991 LLIL 
29 TTHRU=1 
CC TC 54 


11111 CONTINUE 
WRITE(CUT,10CCC) 
WRITE(CUT,10CC3) 
WRITE(CUT,10CC1) 
WRITE{CUT,10C02) 
WRITE(CUT,1C0CC1) 
WRITE(CUT,1C003) 

SSSSS CALL EXIT’ 

END 


1 


22 


12 


WW 
Vad 


aaa 


22 
13 


SLPERECUTINE ANCRMAL(FX,YNCR) 

CIMENSTON X(4C),Y(40) 

COPPEN TE 

CATA X79 Ce .0CO19,G.0062,0.0728,0.C46E8,0.0587,0.2119, 0.350852, 3821 
Ligte4 C2 29 © Sg ODIO ge LTS 9) OI VS Ve TCE 1s OSS Ge OBS Zeke t C2, 
26 S58 Oe SSR Fa Le CGC gL IACI y VI—-3 049-30 9-2 og 2olig—LeDy—l]elig—N ed 
Sir er ge Oia ty = Cw digi o Cig od pl le Sp tin th 9c Big lool gl co Nige will pile Dad elles topple ie Lat 
4SyNPTSPLVS 

PUPP AASAUCK  2ERR CR Sree RPANECS NUMBER, WAS OGUISCE Teri irs te TH 
LENO PUN CIITCh 721K VALE TS) SEY FG Oe v0 8/74 

Be e22 (P=ligNPT 

YE APRS lig L222 

VANCR=Y(T) 

RETURN 

J=Tel 
VACREV(TIF(CYCJI-YOCD) SOX 0S) -X 00) eC FX-XC0T)) 

RETUFN 

CES VESTS 

WRTTECTO,1) 

VN@R=S oN) 

RETURN 

ENG 


SUPROCLTINE GETVAL (NPT XN ,AMEAN,SIGMA,ASJMEAN,SIGAJ) 
VINESR ENTEPRPELATICN FS USEE TE CATAIN VALUES. 
LIMENSTON X02) 9Y02) 9702) W102) gWw(2) 

CCOVYON TC,FN( 200) 

PERVATUY / LC eRe eR=————— VALUES NCT CBTAINES IN SUBRUUTINE GETVAL.§ 
L/2E%, *CLATS SPECTFIED WAS. CUISITILTE TRE FUNCTION LIV LITS.."7 7) 
OE 22 Te NP T,5 

IF (¥N-FACT))11,12,22 

AME AN=FN([41) 

SIGMA=FN(T4+2) 

AJMNEAN=FN(T43) 

SIGAJ=FN(1I+4) 

RETURN 

Pes Gl lS 

IF (FN(T4LI-FN(C14+6)7014512,14 

PCIK UP CCRRECT SUPSRCIPT 

Ee 33a Jels2 

K={-—(S#{?-J)) 

KK=K+4] 

KKK=K+4+2 

K4=K43 

K5=K+4 

X(J)=FNCUK) 

Y(J)=FA (KK) 

Z(J)=FNCKKK) 

Ww(J)=FN(K4) 

WHC J)=FN(KS) 

FACTER=TXN-XOLIIS(IX02)-X(1)) 


NChk MAKE IATERPCLATIOCN 


AMEAR=Y¥(1)4(Y(2)-Y¥C1)) #FACTOR 
S1GMA=Z01)4+(Z(2)-Z(1)) #F ACTOR 
AJMEAN=W(2) 

STCAJ=hW(2) 

RETURN 

CCATINUE 

¥RITECINO,1L) 

RETURN 

EAC 
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SUBROUTINE GRAPH(N,FREC,TITLE) 


SUBRFUTINE TG PLOT CAR GRAPHS 
MAX. CF TEN GARS AVAILABLE 


INTEGER CLT,FLANK ASK 

CIMENSTION FREC(N) »TITLE(2C),ISFREOC( 10,3) sLINE( 88) »NSW(10) 
CATA BLANK,ASK,LINEZ? t,t at ,eRat—'/ 
FCRMAT(IR1L»]CXe*FREGUENCY PLOT ENTITLEN ",290A4/10X,_,"*CANNGT GLE PLOT 
LTE AS NEF. TFOFREQUERTY INTERVALS RECUESTREDD EXCEEDS: PIMitv OFMeROGR 
PAM .*) 

FERMATCULET »L EX, *FREGUENCY PLEOCT TITLE — ".2¢6A447) 
FCRMAT(1CXy,88AL) 
FORMAT(1L2X 914 —92Xe1C(5X%_3N1)) 

FERMAT CLIX, "EACH ©8s Als" “EQUALS ©, 25 "PO TNDSae 
FORMAT(LIX,*CLASS®) 
FORMATIVIX,* INTERVAL *, 1Q0(5X%,12,1X)) 
ECRMAT CEEX, *FREGUENCY * ,10( FS, 3%) 

OO ae 

Die fad beg | eee i) 

KRITECLUT,LITITLE 

GC TC 939SS9 

SEARPGH ECR PAX. 

EPAX=0..5 

ioe Si lig. N 

ASW (T)=9 

JEREC UPS? M=EREE(I ET) 

If CEPREGCT)—EMAX) 22922513 

FRAX=FREE (TI ) 

COATINUEL 

WRITE(CUT, 20 T1TLE 

SET UP FER SCALING IF NECSSARY 

ISCAL=I 

IF (FMAX-50.9)14,14,15 
ISCAL=(FMAX449.3)/59.9 

WRITE (CUT,YS) ASK, ESCAL 

NUR=E+4+R EN 

KRIFE(CUT,.&8) (JFREG(I,1),I=1,N) 
WRITFE(CUT,2).(LINECT).1=1,NUPM) 

CLEAR ARRAY 

UG tls bas 

Ce Vg UA s3 

JFREC( IJ, 1) =RLANK 

SCAL=IESCAL 

MAX=FMAX/SCAL 

fC 44 I[=1,MAX 

X=MAX—(I-1) 

PL 55) =P sh 

LEP UPREOHD) 7 Son l—=X 556 164116 

EGS G6) k= i's 

JFREC(I2K)=ASK 

LE ONS WER) Vile Os Lt 

JFREC(I,2)=PLANK 

NSW(J)=J 

CEOATINUES 

TX=XeSCAL 

WRITECCUT 4) 1X, (( JFREQ( J eK) »K=153),J=1,N) 
fC 88 J=l,yN 


JFREL (Je l= 

WRITS (CUT, 2) (LENE (J) pJ=1_ NUM) 
WRETE(CUT,7) (SEREG (S91) 9 J=1yN) 
WRITE(CUT SE) 

RETUIN 

EAP 


rR 


An 


SLERCUTTNE ZERCLT(N,FREG) 
CIMENSTON FREC(N) 

PERG CVA eer REGUENC LES 
CC 22 T=1,h 

FREC(LT)=C.C 

RETURN 

FAC 


SLPRCUTIANE RAACUCIX,IY,YFL) 
TY¥=1%#6553S 

ITF (1TV)5 9696 
TY=1¥4#721474°3€4741 

YFL=TY 

YFL=YFL* .4E€5€613E-9 

RETURN 

EAC 
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